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Abstract. We investigate an example of noise-induced stabilization in the 
plane that was also considered in (Gawedzki, Herzog, Wehr 2010) and (Birrell, 
Herzog, Wehr 2011). We show that despite the deterministic system not being 
globally stable, the addition of additive noise in the vertical direction leads 
to a unique invariant probability measure to which the system converges at 
a uniform, exponential rate. These facts are established primarily through 
the construction of a Lyapunov function which we generate as the solution to 
a sequence of Poisson equations. Unlike a number of other works, however, 
our Lyapunov function is constructed in a systematic way, and we present a 
meta-algorithm we hope will be applicable to other problems. We conclude 
by proving positivity properties of the transition density by using Malliavin 
calculus via some unusually explicit calculations. 

1. Introduction 

Stabilization by noise is a mathematically intriguing phenomenon. For instance, 
in the classic example of the inverted pendulum, the addition of noise opens up 
a small neighborhood of local stability around a deterministically unstable fixed 
point [AK83 , Git05 ; in the striking examples of [Sch93] , the addition of noise leads 
to global stabilization. In general, however, there are few rigorous proofs of this 
phenomenon for specific systems, and most existing proofs depend upon correctly 
"guessing" a Lyapunov function and then verifying that it satisfies the requisite 
properties. 

In three recent, interesting works [GHW10] IBHWlTj IBD12] . a global Lyapunov 
function is constructed by patching together functions which are locally Lyapunov 
in a collection of regions whose union covers all of the possible routes to infinity. 
These papers are concerned with specific examples in which the stabilization by 
noise is a property of the global dynamics rather than the local dynamics near a 
fixed point. As such, they are closer in spirit to Sch93] than to examples such as 
the inverted pendulum. In these examples, the noise is only important in local- 
ized regions of phase space, but its effect is global, in that it changes the global 
nature of the flow. This nature is hinted at in the patchwork constructions used 
in |GHW10( IBHWIII IBD12] . In |BD12) . this structure is the most explicit; there, 
local asymptotic expansions are used to construct a patchwork of local Lyapunov 
functions. Still in all three works, the local constructions have mainly the flavor of 
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"guess-and-check" with some information of the presumed overall structure of the 
transport in phase space. 

Here we take a more systematic approach to proving global stabilization by noise; 
we outline a meta- algorithm which we hope can be used to produce Lyapunov 
functions in a number of different dynamical systems. Inspired by the examples in 

GHW10, BHW11, BD12], we apply our meta-algorithm to a system of stochastic 
differential equations in the plane whose underlying deterministic dynamics display 
finite-time blow-up for certain initial conditions, but in which the addition of an 
arbitrarily small amount of noise leads to an invariant probability measure. We 
consider essentially the same system as in (ill\\ H)| and one of the examples in 

BHW11 with the specific choice of parameters «i = «2 = 1 and the change of 
coordinates induced by (x,y) > (—x,—y). The choice of parameters is represen- 
tative of the stable regime we are interested in. We demonstrate the existence of 
an invariant measure by constructing a Lyapunov function. Our general approach 
is to build local Lyapunov functions as solutions to associated partial differential 
equations (PDEs), where the PDEs are defined in regions delineated by different 
asymptotic behaviors of the flow. While it is related to that in [GHW10, IBHWll] . 
the Lyapunov function we construct might better be called a "super" Lyapunov 
function, in that it enables us to prove a stronger form of convergence than in 

GHW10, BHW11 . Our exponential convergence results apply equally well to the 
case of degenerate stochastic forcing while those in |BHWll] only prove exponen- 
tial convergence in the uniformly elliptic setting. Our analysis also adapts to the 
specifics of the problem and is likely to produce closer-to-optimal results. 

In DW94| IHKM02| IFMMP081 |Mey08| , the scaling limit of a discrete time 
Markov chain, called the "fluid limit" in the context of these works, is used to 
build a Lyapunov function. In some ways this is related in spirit to our work in 
this paper. However, we are explicitly interested in the case where the noise and 
its fluctuations are fundamentally important to the behavior at infinity. The naive 
fluid limit model, however, is a singular limit and hence does not capture the be- 
havior of those systems in which noise plays an essential role. We will see in our 
example that the naive fluid limit — which is the underlying unperturbed dynamical 
system — is in fact unstable. Our constructions capture the essential stochasticity 
in the regions where it matters at infinity. 

Combining our Lyapunov function with a result on the positivity of transition 
densities, we prove a strong result on the convergence to equilibrium of our specific 
dynamical system. Though the positivity result is neither the most general nor the 
most powerful, it is nevertheless of independent interest, since the proof employs 
sophisticated ideas from control theory and Malliavin calculus in a very concrete 
and transparent way. We hope that it will help develop the readers intuition in 
such matters. 

2. Lyapunov Functions 

We are interested in the stability of a Markov process (X t , Y t ) which is the solu- 
tion to a stochastic differential equation (SDE) on the state space M 2 with generator 
L. In the deterministic setting, a Lyapunov function is a positive function of the 
state space which decreases, often exponentially, along trajectories. In the stochas- 
tic setting, one requires that the function decrease on average. More precisely, we 
define a Lyapunov function V on an unbounded set 7Z C K 2 as follows: 
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Definition 1. A C 2 function V : TZ — > (0, oo) is a Lyapunov function on TZ if 

(i) V{x,y) — > oo as \(x,y)\ oo with (x,y) S 7?., 

(ii) there exist constants m, b > and 7 > such that for all (x, y) £ TZ, 

(LV){x,y) < -mV' r (x,y) + b. 

We say that V is a super Lyapunov function onlZ if 7 > 1 and a STANDARD 
Lyapunov function onTZ i/7 = 1. We call 7 the Exponent of the Lyapunov 
function. If TZ is a strict subset ofM. 2 , we say that V is a local (super /standard) 
Lyapunov function; if TZ is all of M 2 , we say that V is a global (super /standard) 
Lyapunov function. 

We remark that there are several different notions of a Lyapunov function in the 
literature, but the one above will be used in this particular paper. See for example 
}Has80l IMTM1 ISch93l IVer97l IVer99] . 

Remark 2.1. Notice that the continuity of V coupled with its growth at infinity 
implies that the sub-level sets {(x,y) £ TZ : V{x,y) < R} are compact for all R. In 
a certain sense, this is the more fundamental condition, but we will not belabor this 
point here. See, for example, Proposition 5.1 in [HM09 for more details. 

It is well-known that the existence of a global Lyapunov function satisfying the 
properties in Definition [T] implies the existence of an invariant measure [Has80 , 
ISch93[ IMT93| IHM09) . If one adds a mild mixing/minorization condition and as- 
sumes that the Lyapunov function is a standard Lyapunov function, it is possible 
to prove exponential convergence to this invariant measure |MT931 IHM08j . If the 
Lyapunov function is weaker (7 < 1), then the convergence is generally slower 
|MT93| IVer97[ IVer99[ IDFG091 IHaTTO] . in this case, one might rightly call V a sub 
Lyapunov function. 

For the system considered in this paper, we show the existence of a global su- 
per Lyapunov function, along with needed mild mixing/minorization conditions. 
Together, these imply that the rate of convergence to equilibrium is not only expo- 
nential, but also independent of initial condition. 

2.1. General Construction Strategy. We now give an outline of our general 
approach to constructing a Lyapunov function. Since many of the details of the 
implementation depend on the specific example under consideration, this outline is 
meant as an overall rubric. On first reading, this section may seem rather heuristic 
and overly vague. We encourage the reader to take it as motivation at first and 
then reread this section after Section [8] and Section [9] in these later sections, the 
following abstract discussion is made concrete. 

In our general construction algorithm, we begin by identifying a region in phase 
space where there is an obvious choice of a Lyapunov function. We refer to this 
region as the "priming" region; it is often characterized as a subset of phase space in 
which the deterministic flow is directed toward the origin. We refer to the associated 
local Lyapunov function in this region as the "priming" Lyapunov function. Next, 
by contrast, we identify the region in which the deterministic dynamics exhibit 
instability — for example, blow-up in finite time — and for which noise is essential to 
the stabilization, at least insofar as the noise ensures that the system leaves this 
region. Since this region is noise-dominated to some degree, we refer to this region 
as the "diffusive" region. The construction of a local Lyapunov function in this 
diffusive region is a key component of our methodology. We then "propagate" the 
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priming Lyapunov function, from the priming region to the diffusive region, through 
a series of intermediate regions of phase space, which we call "transport" regions, 
until we have covered all possible routes to infinity. Following this prescription, 
we obtain a sequence of local Lyapunov functions, which we mollify to obtain a 
global Lyapunov function. Beyond this general, overarching strategy is a philosophy 
of determining the relevant scaling at infinity in each of the above regions and 
systematically producing local Lyapunov function which respect this scaling. 

To determine more precisely the boundaries of these regions, we study the scaling 
of the generator L of the SDE as \(x, y)\ — > oo. Each region corresponds to a 
different dominant balance of the terms in the generator. (Sec [B0781 IWhilQ] 
for a discussion of the concept of dominate balances in asymptotic analysis and 
Section E] for the details in our setting.) In order to facilitate the mollification of 
the local Lyapunov functions, we choose the regions so that the intersection between 
adjoining regions is both nonempty and unbounded. Neglecting all but the terms 
involved in the associated dominant balance, each region also has a differential 
operator associated to it which captures the dominant behavior of the generator 
in the region as |(a;, j/)| — > oo. Beginning with the region adjacent to the priming 
region, we propagate the priming Lyapunov function through the adjacent region 
by solving an associated Poisson equation of the form 

(Lv)(x,y) = -f(x,y) 

v(x,y) = g(x,y) on the boundary. 

The differential operator L is governed by the dominant balance in the region under 
consideration. Since it represents a dominate balanced at infinity, it is necessarily 
an operator which scales homogeneously. Hence if the righthand side and the 
boundary conditions are chosen to scale homogeneously at infinity in a compatible 
way the solution v will also scale homogeneously at infinity. The boundary data g 
for the Poisson equation is given by the dominant behavior/scaling of the priming 
Lyapunov function on the boundary between the priming and adjacent region. The 
right-hand side, /, of the Poisson equation is chosen to be a positive definite function 
which grows unboundedly and satisfies certain scaling properties that we specify in 
Section 19.31 and that are compatible with the scaling of the region. 

We iterate this procedure to construct local Lyapunov functions as solutions to 
associated Poisson equations in each of the transport regions. Furthermore, we 
construct a Lyapunov function in the diffusive region by solving a Poisson equation 
as well, again with the boundary data determined by the dominant behavior of the 
local Lyapunov function in the adjacent transport region. An advantage of this 
approach, therefore, is its consistency: the same procedure is used to construct 
local Lyapunov functions in all but the priming region (where the local Lyapunov 
function is usually straightforward to deduce) . 

As we solve the sequence of Poisson equations, we encounter boundaries with- 
out boundary data. While a priori this could be an issue, we will see that in 
our model problem, in all but the diffusive region, the deterministic dynamics are 
dominant. Hence the associated Poisson equations are governed by first-order op- 
erators requiring only one boundary condition. This is consistent with the idea of 
the priming Lyapunov function being propagated through a sequence of transport 
regions. Again, a priori this could lead to an incompatibility between two different 
boundaries of a given region, particularly if the relevant operator in a region is only 



PROPAGATING LYAPUNOV FUNCTIONS 



5 



first-order and cannot accept generic initial data on two boundaries. However, in 
our model problem and all of the other problems we have explored, sequences of 
compatible transport regions are separated from each other by diffusive regions. 
Since the associated differential operator in the diffusive region is second-order, 
the associated Poisson equation produces a smooth solution even with all of its 
boundary data specified. 

3. The Model Problem 

As our model problem, we consider essentially the same problem as in (GHWIOl 
BHWlTj and which was suggested to us by one of the authors: 

dX t = {X 2 - Y t 2 )dt + V2o~ dW t {1) 

dY t = 2X t Y t dt + ^/2o~ dW t {2) 

with a x > and o~ y > 0. Notice that when 

{<j Xl /jy) — (0,0), the resulting deterministic equation blows up in finite time if 
xq > and yo = 0. In light of this, it is striking that for any a y > 0, system ((T|) has 
a unique invariant probability measure it. This was first proven in [GHW10 and is 
also a consequence of one of our main results, which is given below in Theorem^ In 
the sequel to [BHWll] . the authors prove exponential convergence to equilibrium. 
The principal difficulty in both of these works was the establishment of a standard 
Lyapunov function. 

Let Pt be the Markov semi-group associated to the process (Xt,Yt) and defined 

by 

(2) (Ptflfav) = U(x, v )[<KXt,Y t )]. 

Define the action of Pt on a probability measure fj, by (fj,Pt)(A) — J P t (x, A)fi(dx) 
for any measurable set A. An invariant probability measure fj, is any measure such 
that iiP t = /i for all t. 

We prove the following theorem, which is stronger than the previously cited 
results on the existence of a standard Lyapunov function. 

Theorem 1. There exists a C 2 function V : M 2 — > (0, oo) which is a super Lya- 
punov function for the dynamics given by {T]). More exactly, for any choice of 
8 € (0, |) the super Lyapunov function V can be chosen to have an exponent §§i§ 
and satisfy c\(x,y)\ s < V(x,y) < C\(x,y)\^ s+ i for some positive constants c and 
C. 

The existence of an invariant measure fj, is an easy consequence of this theorem. 
To determine rates of convergence to the equilibrium measure /i, we introduce the 
following family of weighted total variation metrics. For j3 > and probability 
measures /ii and /i2, we define 

(3) M/ii,M2)= su P / ^( z )(Mi - f-2)(dz) 

\\4>\\0<^J 

where 

Notice that po is just the standard total- variation norm. 
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The standard Lyapunov function and supporting estimates developed in [BHW11 
essentially show that there exists positive C and ij so that 

pi(jjLiP t , /i 2 -Pt) < Ce-^pxim,^) 

for any probability measures p\ and pi. Using Theorem [T] on the existence of a 
super Lyapunov function, we establish the following stronger convergence result: 

Theorem 2. If a y > 0, then for any (3 > there exist positive C and r\ so that for 
all probability measures p\ and pi one has 

pp(ji,lP t , fjt 2 Pt) < Ce-^Wfii - p 2 \\rv 

for all t > 0, where \ \tv represents the total variation norm. Here the constant C 
depends on the choice of ft but the constant r\ does not. 

The strength of this result is that the dominating norm on the right-hand side 
is scale and translation invariant. As we will see, when a super Lyapunov function 
exists, one can usually prove a stronger result than the standard Harris- type ergodic 
theorem associated to a standard Lyapunov function. 

Remark 3.1. As already mentioned, the existence of an invariant measure p fol- 
lows quickly from Theorem^ or Theorem^ The fact that there is only one invariant 
measure is immediate from Theorem \^ 

One consequence of our estimates is the following information on the unique 
invariant measure p. The proof of Theorem [3] below is given in Section 111. 31 

Theorem 3. As long as a y > 0, then p has a smooth density with respect to 
Lebesgue measure which we denote by m(z). If o~ Xl a y > 0, then m(z) > for all 
z G M 2 . If a x — and a y > 0, then m(z) = if z = (x, y) with x > 0, and m(z) > 
if z = (x, y) with x < 0. 

4. Outline of Paper 

In Section |6l we show how the existence of a super Lyapunov function leads to a 
strong regularization of moments. In Section [3 we discuss further the properties of 
the deterministic model problem. In Section [51 we perform an asymptotic analysis 
of the generator associated with ([1]). In Section we use associated Poisson equa- 
tions to construct local super Lyapunov functions in the different regions whose 
boundaries are determined by the asymptotic analysis. In Section [TUl we patch the 
local Lyapunov functions together to construct the global Lyapunov function and 
thereby prove Theorem Q] In Section HT1 we prove, under various assumptions, the 
existence of a smooth transition density with various positivity properties. Our ap- 
proach here invokes methods from geometric control theory and Malliavin calculus 
in a manner which, we hope, will be of independent interest. In Section lll.3[ we 
transfer the smoothness and positivity results to the invariant measure and in doing 
so prove Theorem [3] In Section [T21 we prove that Theorem [I] when combined with 
a standard minorization condition, implies Theorem [5J In Section 112. 1[ we show 
how in the uniformly elliptic setting, namely a x ,a y > 0, the needed minorization 
condition follows immediately from the strong from of positivity which holds in 
that setting. In Section 112.21 we show how the weaker positivity properties which 
hold when a x = 0, a y > are sufficient to prove the minorization condition. In 
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Section I13[ we make a few concluding remarks. The Appendix contains a rela- 
tively standard comparison result which we include for completeness. It is used in 
Section |B] about the Super Lyapunov Structure. 
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6. Consequences of Super Lyapunov Structure 

We begin with a lemma, whose proof is given at the end of the section, which is a 
simple translation of the bound on the generator for the definition of a global super 
Lyapunov function to a bound on the action of the semigroup. Despite its simplicity, 
it is nonetheless the key to all of the enhanced results that are a consequence of the 
existence of a super Lyapunov function (as opposed to merely a standard Lyapunov 
function) . 

Lemma 6.1. Suppose that V : R 2 —¥ (0, oo) is a super Lyapunov function for the 
SDE corresponding to a Markov semi-group P t . Then for every t > 0, there exists 
a positive constant K t , such that t i— >■ K t is a continuous, monotone decreasing 
function on (0, oo) with K t — ¥ (26/to) 1 / 7 as t —¥ oo, and 

{P t V)(z) < K t for all zeR 2 and t > . 

Recalling the definition of pp from @, Lemma 16 . 1 1 implies the following result. 

Proposition 6.2. IfV is a super Lyapunov function and K t is the constant defined 
in Lemma \6.1\ then for any t > 0, /3 > 0, test function <f>, and probability measures 
pi and fi2, we have 

\\PtHo < (1+ PKJWWp and pp(p 1 P t ,p 2 P t )<(l + f3K t )\\p 1 -p 2 \\ T v- 

Remark 6.3. It is clear that po(Ml> M2) — \\Pi — A^IItu and furthermore if < a, 
(3 > and K = sup^ iljl^y^j , then one has \\4>\\p < -fC||i?!)|| a , which implies 

{<t> ■■ H\\ a <l/K}c{<j>: UWp <1}, 

which in turn implies p a (/^i , ^2) < Kpp^pi, (12)- Thus as long as a and (3 are both 
positive, the associated norms and metrics are equivalent. However, if one of them 
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is zero, the needed inequalities only go in one direction. Nonetheless, Proposition ^. 2\ 
allows us to use the action of Pt to recover the missing inequality. 



Proof of Proposition \6.2\ By similar reasoning to that used in the second part of 
Remark 16.31 we see that if one assumes that ||P t ^|| < (1 + /^-K^H^H^ for some 
constant K u then j> : \\<j)\\p < 1/(1 + f3K t )} C {<j> : \\P t <j>\\o < 1} which then 
implies that pp(piP t , ^Pt) < (1 + PKi)po(^i, P2)- Since as noted in Remark IB~3l 
Po(pi, A*2) = 1 1 /-ii — 1| j the proof of the second quoted inequality is now complete 
provided we prove the first. 

Now since \<f>(z)\ < \\<f>\\p{l + (3V(z)) for all z, one has 

\(Pt4>)(z)\ < Uh(l + f3(P t V)(z)) < 1101^(1 + ^). 

Since the right-hand side is independent of z, we obtain the desired result by taking 
the supremum over z. □ 

Remark 6.4. In light of Remark \6.3\ and Proposition 1 6. 2\ to prove Theorem^ we 
need only prove the more standard Harris chain-type geometric convergence result 
of Pp{niPtiP2Pt) < Cexp(-r?t)p^(/ii,/x 2 ) for some (3 > 0. 

Proof of Lemma \6.1\ Let V t — V(Z t ), where Z t is the solution to the SDE corre- 
sponding to Pt- Let L denote the generator associated to the SDE corresponding 
to Pt- Since V is a super Lyapunov function, there exist constants m, b > and 
7 > 1 such that 

LV t < -mV t y +b for all t > . 

By Dynkin's formula, 



(P t V)(z) = E z [V t ] = V(z)+E z 



LV.ds 



ii 



< V(z) - m 



/ E,[V?]da + bt 
Jo 



< V(z) - m / E^Vs} 1 + bt by convexity. 



For simplicity of notation, let <fi z (t) — (PtV)(z) — E z [Vt]. Then 4> z (t) satisfies the 
following differential inequality: 

<t>' z (t)<-m[ct> z (t)V + b 

1 

<-j[Mt)F if M*)>Q 

Let R = and let r = inf{t > : <j) z {t) < R}. Since <j)' z (t) < if <p z (t) > R, 

this implies that once (f> z (t) < R, 4> z (t) remains less than or equal to R for all times 
afterward. Thus, for all t > r, 4> z (t) < R. Now suppose ip z (t) satisfies the following 
differential equation: 

777 

Then by Proposition IA. II in the Appendix, (f> z (t) < ip z (t) for all t € [0,r]. Now the 



differential equation for ip(t) can be solved explicitly to obtain that for all t € [0, r] 
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Defining the constants Kt as follows 




we conclude that 4> z (t) < K t for all t > 0, which completes the proof. □ 

7. Deterministic Equation 

To better understand the context of our results for the stochastically perturbed 
system, we pause for a moment and highlight some properties of the underlying 
deterministic dynamics: 

(4 ) it = ^t-Vt 

y t = 2x t y t ■ 

The trajectories of the system are shown in Figure [TJ from which the dynamics of 
the system can be quickly and easily understood. 

For any initial condition (xo,yo), the solution (x t ,yt) to this system is given by 

~ {xj + Vp)t 

(5) Xt (i - x ty + ( yo ty 

yt 



(l - x t) 2 + (yot) 2 ■ 

In particular, the system exhibits finite-time blow-up (at time t = — ) for initial 
conditions (xq, 0) on the positive x-axis. For all other initial conditions, the w-limit 
set w(xo,j/o) is simply the origin, which is the unique fixed point of the system. 
We note that the origin is not reached in finite time by any trajectory with initial 
condition (x ,y ) ^ (0,0). 

Now, for any choice of initial condition (xo, yo) not on the rr-axis, the trajectories 
of the deterministic system are circles centered at the point C(xo,yo) with radius 
R(xa,yo) given as follows: 

(6) C(xo,yo) = (0,4±% R(x ,y ) = ^f. 

v 2y / 2\yo\ 

Furthermore, for all choices of initial conditions (xo,yo) not on the positive x- 
axis, the time to return to a fixed ball of radius R about the origin is uniformly 
bounded by j| . In Section 112.21 we employ this uniform bound to prove a posi- 
tivity and minorization condition on the transition density for the stochastically- 
perturbed system. 



8. Dominant Balances of Generator 

We now begin the program laid out in Section [2. II We begin by considering the 
dominant operators in various regions of the state space. 
Associated to the SDE ([TJ is the generator L defined by 

(7) L = (x 2 - y 2 )d x + 2xyd y + a x d xx + cr y d yy . 

In order to prove that the addition of noise arrests the blow-up on the x-axis 
sufficiently to produce an invariant probability measure, we need to understand the 
behavior of the dynamics at infinity. There are many different routes to infinity and 
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Figure 1 . A number of representative orbits of the deterministic 
dynamics governed by ((4]). 

we now consider the various possible dominant balances associated with different 
routes. 

To help identify the relevant scaling, consider the behavior of L under the scaling 
map (x, y) i— > (£x, £ p y) which produces 

£ x 2 d x - £ 2p ~ l y 2 8 x + £ 2xyd y + £~ 2 a x d xx + r 2p a y d yy . 

If p = 1 the first three terms balance and dominate the remaining terms as £ — > oo . 
If p > 1 then the second term dominates. If p = — \ then the first, third and 
fifth balance and dominate all other terms as £ — > oo. These balances cover all 
of the routes to infinity except for those which approach or rest on the y-axis and 
identify p = —1/2 as a critical scaling. (The routes near the y-axis are captured by 
p = —1/2 and £ — > but these will not play an important role in our analysis.) 
If \x\y 2 < oo as \(x, y)\ — > oo with x > a > 0, the dominant part of L is contained 

in 

(8) A = x 2 d x + 2xyd y + a y d yy . 

If \x\y 2 — > as \(x, y)\ — > oo with x > a > 0, then the dominant part is only d yy . 
Notice that d yy is contained in A, so we can still choose to use A this region. In all 
other relevant cases as \{x, y) \ — > oo, the dominant part of L is contained in 

T = (x 2 - y 2 )d x + 2xyd y . 

We have neglected the term a x d xx in the operator T which scaling analysis suggests 
might be relevant in neighborhood of the y-axis. However its inclusion does not 
qualitatively change the behavior in a neighborhood of the y-axis. The same can 
not be said of the term <J y d yy in a neighborhood of the ir-axis. 
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8.1. Scaling. To better understand the structure of the solutions in the various 
regimes, we investigate the scaling properties of the various operators introduced 
in the previous section. We introduce the scaling transformations 

Sf> : (x, y) ^ (tx, riy) and sf ] : (x, y) ^ (tx, ly) . 

Observe that operator A scales homogeneously under the scaling Si , while the 

(2) 

operator T scales homogeneously under the scaling SI . We would also like the 
operator T to scale homogeneously under the scaling Si^; however, this does not 
hold for all of the terms in T. We remedy this by introducing a non-negative 
parameter A and defining the family of operators 

(9) T\ = [x 2 - Xy 2 )d x + 2xyd y 
and extending the definition of the scaling operators by 

Sjp: (x,y,X) H- (tx,£-iy,t 3 X) and S^p : (x, y, A) i-> (tx, ly, A) . 

Now T\ scales homogeneously under the scaling map Sg and A remains invariant 
under . This gambit of introducing an extra parameter to produce a homoge- 
neous scaling was also used in a similar way in [CMMSlTj . 

Given a function cf> : 1Z x [0,oo) — > M, where 1Z C K 2 , we say that (j> scales 
homogeneously under the scaling Sf 1 ' if <f> o = £ 5 (f> for some S. In this case, we 
say that <j> scales like £ s under the i-th scaling. We write this compactly as <f> ~ £ s . 

Proposition 8.1. If <f> ~ I s then d x (j) ~ £ & ~ x and d y <j> ~ £ s+ i . Similarly, if (j) ~ £ s 

then d x <j) ~ £ S ~ 1 and d y <f> ~ £ s ~ 1 . In both cases, if one side is infinite, then so is 
the other. 

Proof of Proposition \8.1l We only show one case; all others follow similarly. If 
(j) ~ £ 5 , then cf>(£x,£~iy,£ 3 X) — £ 5 4>(x,y, A). Differentiating in x, we obtain 

l(d x $)(tx,t-^y, ex) = £ 5 (d x cf)(x, y, X) . 
Dividing through by I, we conclude that d x <p ~ £ s ~ 1 . □ 

In the next section, we decompose the plane into regions where the different 
dominant balances hold. These regions are defined by boundary curves which are 
well-behaved under one or both of the scalings. To facilitate the construction of 
these regions, given xo > 0, yo > 0, A > and p € R, we define the following 
"elementary" regions: 

a / u f + Xy 2 ^ xl + Xyl \ 
K(x , y , A) = \ n > \ 

(10) r±(as , W>) = {±x > xo, \x\ p \y\ < x p y } . 
Observe that for any £ > 0, we have to following scaling relations 

sV(T±(x ,yo))=T±(tx ,t-ho), S ( / ] (T±(x Q ,y )) = r±(£x ,ty ), 
SP(A(x ,y ,X)) = A(tx Q ,£-iy ,£ 3 X), S*p (A(x , y , A)) = A(£x , £y , A), 
and lastly Tf (£x , £~ p y ) C r±(a;o,yo) for £ > 1 . 
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X 



Figure 2. The different regions in which local Lyapunov func- 
tions are constructed. 1Z\ is the priming region. The two regions 
labeled IZ2 are transport regions. And IZ3 is the diffusive region 
which connects the two transport regions in which information is 
propagating in different directions. 

9. Construction of Local Lyapunov Functions 

Based on the discussion in the previous section, we will divide the plane into 
three regions lZi(a), where a is a positive parameter that we specify later. As 
described in Section 12.11 we call these regions the "priming," "transport," and 
"diffusive" regions, respectively. We now describe the placement of these various 
regions which are indicated pictorially in Figure [2j 

Our priming region, 1Z\(a), is a subset of the left-half plane, and here there exists 
a very natural Lyapunov function, because in this region, the deterministic drift 
is directed toward the origin. On the other hand, the diffusive region, TZs(a), is a 
funnel-like region around the positive x-axis where there is finite-time blow-up in 
the deterministic setting. Demonstrating the existence of a local Lyapunov function 
in the diffusive region is a key piece in proving noise-induced stabilization in our 
model problem. The transport region IZ2 (ot) is governed primarily by deterministic 
transport from the diffusive region to the priming region. In this section, we focus 
on the construction of a local Lyapunov function in each of these three regions. 

9.1. The Priming Region. When looking for a priming Lyapunov function, it is 
natural to consider the norm to some power. In this specific example, we expect 
the norm to some power to be a Lyapunov function in the left-half plane since the 
drift vector field points at least partially towards the origin; see Figure Q] 
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For 8 > 0, we define V\{x, y) — (x 2 + y 2 ) 2 and observe that 



L Vl (x,y) =8x(x 2 + y 2 )U 8(°- - l)(x 2 + y 2 )^ 2 (2a x x 2 
(11) i 

+ (cr x + a y )5(x 2 +y 2 )?- 1 . 



In particular, if (x,y) £ r_ 1 (#, 1), we get that 
a 8 



2<r y y 2 



{Lv!)(x,y) < - 



\Ju 2 + 4 



(x 2 +y 2 )^ + 8(8-2)^ 



a x x 2 + o y y 2 



(x 2 + y 2 ) 



2\2- 



(x 2 + y 2 ) 



2>1- 



q5 



Va 2 + 4 



(^ 2 + z/ 2 ) i 



Va 2 +4 ^[ ^-2] (<t x x 2 +(i y y 2 ) 
(x 2 + y 2 )^ 



(x 2 + y 2 ) I 



This implies that for any 5 > and a > 0, there exists an i?i sufficiently large so 
that if \(x, y)\ > Ri, then the term in the square brackets is greater than i. Hence 
v% is a super Lyapunov function in the region 

^i(a) = r: 1 (|,i) 

with exponent As we will see later, we will have to restrict S to the interval 

(0, |), and this automatically implies that S G (0,2). In turn, this guarantees that 
S — 2 < and that the term in the square brackets above is greater that i provided 



1 

< r- 



a 



(x 2 +y 2 )2 2 ^Ja 2 + 4 
We formalize this observation in the following proposition. 

Proposition 9.1. For any a > and 5 € (0,2), if (x,y) £ lZi(a) with \(x,y)\ > 
Ri , then V\ satisfies 

{Lvi){x,y) < -mi vj^x.y) 



where mi — 



a.6 



> 0, 71 = > 1, Ri = 2 (a x + o- y ) 



Our choice of the region IZi(a) is motivated by the following. From (fTTj) . it is 
clear that we need to define a region in the negative half-plane bounded away from 
the y-axis. Furthermore, in order to guarantee that v\ is super Lyapunov, we need 
to ensure a region in which |(x,y)| — > 00 implies \x\ — > 00. Note that 1Z\{a) is a 

subset of the left half-plane, in which the dominant dynamics at infinity are given 

(2) 

by T and hence the relevant scaling transformation is S e . For this reason, it is 

(2) 

desirable to define the boundary of the region so that it behaves well under . 
From the previous section, we see that 



SP(rz 1 {x ,y )) = rz 1 (£x 0> £y )cTZ 1 (x ,y ) for l>\ 



which motivates our choice of 7?-i(a) and the shape of its boundary in particular. 
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9.2. Decomposition of Remainder of Plane. We will now propagate the prim- 
ing Lyapunov function through a sequence of regions until all of the routes to 
infinity are covered. 

As mentioned above, near the boundary of 72i(a) and away from the x-axis, 
the operator T is dominant. This holds true until one enters the region defined by 
the curves xy 2 = c where c is a sufficiently large positive constant and x > is 
sufficiently large. At this point, the dominant balance changes and the operator 
A becomes dominant. Hence we define the transport region, 722(a), with one 

(21 

boundary inside the region 72i(a) which is invariant under the scaling S\ , and 
one boundary which is defined by the curve \x\y 2 — c for some constant. As we 
make precise in the definition below, we will choose c = a. 
We set 722(a) = 7^2(a, 1) where for a, A > 0, we define 

72 2 (a, A) = rf (ax/A, l) c n A(ax/A, 1, A) n TZ^ay/X, l) c . 

2 

Now, observe that outside of 72i(a) U 722(a) all of the routes to infinity have 
\x\y 2 < 00. Hence the operator A is dominant in this entire region and we do not 
need to further subdivide the remainder of the plane. To define IZ3 (a), recall that 
we need nontrivial overlap with the transport region 722(a). Hence we again chose 
a boundary curve of the form \x\y 2 — c but with c > a. In particular, we define 

72.3(a) = Ft (2a, 1). Note that 72.3(a) is the diffusive region: the diffusion term in 

2 

the operator A is critical to the stabilization of the process here. 
In summary, for each a > 0, we have defined three regions 

Priming Region: 72i(a) = TZ 1 (^, 1) 

Transport Region: 72.2(a) = 722 (a, 1) 

Diffusive Region: 723(a) = Ft (2a, 1) . 

2 

Notice that 72 1 (a) n 722 (a) and 722 (a) fl 723 (a) are nonempty and that R 2 \ (72 1 (a) U 
722(a) U 723(a)) is a bounded set. We construct a local super Lyapunov function 
in each of the three regions and then smoothly patch them together to form one 
global super Lyapunov function on the entire plane. 

9.3. The Associated Poisson Equations. We now propagate the priming Lya- 
punov function v\ which is defined in 72i(a) to the other regions by solving a 
succession of Poisson equations. Throughout most of our construction, a will re- 
main a free parameter; we specify a later to ensure a number of necessary estimates. 
We begin with the transport region 722(a). 

9.3.1. The Transport Region 722(a). For 8 > and a > 0, we define v 2 (x,y) as the 
solution to the following Poisson equation: 

/ 2 1 2\< 5 + 1 

(12) J (Tv2)(x,y) - - (^ L ) ^ 72 2 (a) 

v 2 {x,y) = vi(x,y) on dBi(a) 

where dBi(a) — \ x < —a, lyl = — \x\ \. 

I a J 

The rationale for this is as follows. We wish to propagate the priming Lyapunov 

function through the region 722(a), so we need to take it as the boundary condition. 

Since the operator T represents one of the dominate balances, it necessarily scales 
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homogeneously. In this case, T scales like t 1 under the scaling transformation S\ . 
Hence if i>2 is to scale homogeneously under as £ p for some power p then the 
righthand side must scale like £ p+1 and the boundary conditions must scale like l p 

C2) (21 

both under . (Notice the boundary dB\ is invariant under S\ .) 

Notice that our choice of right-hand side scales as £ s+1 and the boundary con- 

<- (2) 

ditions scale as would be consistent with a solution which scales like I under S\ . 
The form of the boundary conditions are dictated by our choice of v\ . The exact 
from of the righthand side was chosen so that it was constant along trajectories of 
the limiting dynamics in 722(a) which are the characteristics of T. 

9.3.2. The Diffusive Region 72.3(a). For S > and a > 0, we define V3(x,y) by the 
following Poisson equation 



and c\, C2 > are constants which will be chosen later. We remark that the values 
of ci and ci do not affect the local super Lyapunov property of V3, but rather are 
chosen in order to facilitate the patching of the local super Lyapunov functions into 
one global super Lyapunov function in Section fTOl As before, we have chosen a right- 
hand side which is negative definite, scales homogeneously under the appropriate 
scaling, namely , and has unbounded growth in the region. We use a constant 
multiple of x s as the boundary condition rather than the function V2 from the 
neighboring region because we want a function which scales homogeneously under 
. However, x 6 is in fact the asymptotic behavior (up to a constant multiple) of 
V2(x,y) as \(x, y)\ 00 on the specified boundary. 

In Section 19.51 we verify that V2 and V3 are super Lyapunov functions in the re- 
gions in which they are defined. However, we first establish a number of preliminary 
results. 

9.4. Existence of Solutions and Their Properties. The scaling properties of 
the solutions to the above Poisson equations are one of main tools we use to show 
that they are local Lyapunov functions. This is because, with one exception, points 
at infinity in a given region can be scaled back to points in the same region by the 
scaling transformation under which the associated differential operator is homoge- 
neous. As we discuss below, the exception is the subregion of 72.2(a) which lies near 
the boundary of 72.3(a). 

9.4.1. Properties of the Solution in the Transport Region. Care must be taken when 
scaling the points in the subregion of 722(a) which lie close to the boundary of 
723(a). The points in this region naturally scale with while the operator T 

(21 

which is associated to 722 scales homogeneously under . This issue was also 
addressed in Section [8.11 where we introduced the parameter A to generate a family 
of operators T\ which scale homogeneously with . 




on 723(a) 
on 9i?2(a) 



(14) 
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With this mind, it is natural to introduce the function V2(x, y, A) which, for 
a given A G (0,1], solves the following family of auxiliary Poisson equations in 

(T\v 2 ){x,y, A) = ~h(x,y, A) on 72 2 (a, A) 
v 2 (x,y, A) = f(x,y, A) on dB^aVX) 




/(z,j/,A) = A— (x 2 + A 2 / 2 )^ 



12/1 

For ease of notation, we write 

h(x,y) = h(x,y,l) and f(x,y) = f(x,y,l). 
Notice that ft ~ / ~ ^, h ~ and / ~ ^ where 5 was defined in (14 



Also observe that V2(x,y, 1) coincides with the v 2 (x,y) defined by ([12]) . 

9.4.2. Properties of the Solution in the Diffusive Region. The dynamics associated 
to the operator A, which is dominant in 723(a), should be understood as having 
one diffusive direction and one deterministic direction which is uncoupled from the 
diffusion and acts as the "clock" of the diffusion. To see this, observe that A is the 
operator associated to the system of SDEs given by 



dX t = Xfdt X 



x 



(16) 

dY t = 2X t Y t dt + y/2a y dW t Y a = y . 

Now, let (X ,%) = (x,y) lie in 72 3 (a) and define f = inf{t > : (X t ,Y t ) € 
dB2(a)}. Then 113(3:, y), which was defined in (|13l) . can be represented probabilis- 
tically as 



"3(2:, y) = C2l(x,3/) [X?] + ciE (Xi 



X^ds 



(17) =(f + c 2 )E {x<y) [X*]-fx s 

provided that, first, the expectation is finite; and second, that the right-hand side 
of equation (fl"7|) depends in a C 2 fashion on (x,y) £ 72.3(a). Both of these facts 
will follow from Proposition 19.21 which we present below, and are made formal in 
Proposition I9.3[ which appears in the next section. 

Since X t is deterministic, this representation of V3 amounts to a deterministic 
function of f. To better understand the properties of f, we introduce the time 

change T(t) — J Q X s ds — — In |1 — xt\ and the process Z T ^ = X^ Y t . Due to the 
scaling of the boundary of 72.3(a), if we define r = inf{T > : \Zt\ > V2a} then 
f = i(l - e~ r ), X t = xe Ti - l \ and Z T satisfies the SDE 

(18) dZ T = ^Z T dT + dW T , Z = x?y. 

Since Zt is the solution to a Gaussian SDE, the following proposition follows easily. 

Proposition 9.2. For 5 < | and (x,y) g 72.3(a), y ^ [e 5r ] < 00 and the map 

{x,y)^E (Xiy) [e* T ] isC 2 . 
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Proof of Proposition \9. e A To see the finiteness of the expectation, observe that 
P {x , y) (e* T > S )=P {X J sup \Z T \<V2aj <F {Xty) (\Zi lns \<V2aj 



V 0<T< lns 



■ In s 



{j2c~s^\ S e-'l r dWr < y/2a) < ( ^ V 

V V v In J ~ \a y K( S s/ s -iy 



Hence for 8 < |, this decays sufficiently rapidly in order to guarantee that E( x>y ) [ e<5r ] 
is finite. The continuity properties now follow from the continuity properties of r. 
Specifically, E^j,) [e 5r ] = g(y/xy) where g(z) solves the following ordinary differ- 
ential equation 



(19) 



<J y g"{z) + \zg'(z) + 8g(z) = for g £ (—/2a, \/2a) 



[g{V2a) = g(-V2a) = 1. 
Since by standard results on the regularity of ODEs, g(z) S C 2 ([~y/2a, via]), we 



conclude that E^^) [e ST ] — g(\/x y) € C 2 (IZ3 (a)) as desired. □ 



We remark that this proposition imposes a further restriction on the size of the 
parameter S, which previously was only required to be positive. Observe that in 
light of (H3J) the requirement that 8 < | forces 8 £ (0, |). 

9.4.3. Principal Result on Existence and Scaling of Solutions. We consolidate these 
observations and now state and prove our principal existence and scaling result. 

Proposition 9.3. For every 8 £ (0, |), there exists a strictly positive C 2 function 
v 3 : ~R-3( a ) ~^ (0,co) which solves (|13p . For every A £ (0,1], there exists a strictly 
positive C 2 function V2 '■ TZ-2 (pi, A) — > (0, 00) which solves (|15j) . In addition, V2 ~ I s , 
v 2 ~ £ s , v 3 ~ £ s and (x,y,X) H> v 2 (x,y, A) is continuous on TZ^(a) x [0,1] where 
lZ 2 (a) = ^\£[o^i]R-2(a, A). In fact, v 2 has an explicit formula given in (1211) below 
and V3 a semi-explicit formula given in (|20p afeo below. 



Proof of Proposition \9.3l We begin with U3 . The preceding discussion all but gives 
the existence proof. In particular, it shows that if g is defined by (|19j) and 8 by 
([14]) then the map 

(x, y) ^ E (Xiy) [xf] = x% Xiy) [e Sr ] = x S g(^y) 

is well-defined, positive, and C 2 for 8 £ (0, |) and (x,y) £ IZsia). Returning to 
(TT7)) . classical results (see, for example, Bas98 ) allow us to justify the stochastic 
representation formula for v% , which now can be rewritten as 

(20) v 3 (x, y) = x S [(f + c 2 )E (x , B) [e* T ] - f ] = J [(f + c 2 )g(^y) - f 

As a consequence of this formula, to prove the scaling it suffices to show that 

E(x )V ) [e Sr ] ~ £° ■ 

This is clear, since the only dependence of r upon x and ?/ results from Zq, and 
= x"Sy = (£x)i (y£~^) is invariant under Sj^. We now turn to v 2 . While in 
light of the scaling and continuity properties of T\, h, and / this result could be 
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obtained by abstract means, we employ the method of characteristics to produce 
an explicit solution. Namely, 

(2i) Bfes , A) = (f!+V)«[ 1 f I+(l ^ + ^(_J_) 

The scaling properties, regularity, and positivity follow by direct calculation with 
(ED). □ 

Remark 9.4. As A — > 0, V2(x, y, A) given in (|2"Tj) converges to x 2S+1 \y\~( s+1 K This 
was expected since formally taking A — > in (|15p produces the equation 

x 2(6+i) 

())) (x 2 d x + 2xyd v )v = ~^pr ™ft 2 (a,0) 

f(a;,y) = on dP>i(aV\) . 

The solution to this simplified equation is easily seen to be x 2S+1 \y\^^ +1 ^ . Even in 
a setting where (|15l) cannot be solved explicitly, this simplified equation may well 
be easier to solve. We will see in Section 1101 that the most delicate parts of the 
patching require precise information about the limit of V2 as A — > 0. This suggests 
that the analysis may be feasible even when (|15p is not solvable. 

Remark 9.5. For both V2 and we have used specifics of the solutions to verify 
the scaling. It is possible to derive the results just from the scaling of the operators, 
right-hand sides, and boundary conditions. The positivity for both solutions also 
follows from the positivity of the boundary data and the negative definiteness of the 
right-hand sides. 

9.5. Proof of Local Super Lyapunov Property. Letting Br(z) — {(x,y) € 
M 2 : \(x, y) — z\ < R}, we state the following proposition which establishes that V2 
and are local super Lyapunov functions. 

Proposition 9.6. Fix any 5 € (0, |) then for all a > sufficiently large, there 
exist constants rrii > and Ri > so that if {x, y) € 7?.,(a) with \ (x, y)\ > Ri, then 
Vi satisfies 

(Lvi)(x,y) < -rriivf {x,y) 

for i = 2,3 where 72 = 73 = ff+f- ^ n addition, we have the following refined 
estimate in the second region which emphasizes its transitional nature and which 
will be of later use. Defining 



(23) 7l 2 1] (a) = 112(a) (a) c and (a) = K 2 (a) D T+^a, l) c , 



if j = 1,2 and (x,y) € T&P 1 we have 



(Lv 2 ){x,y) < -m 2 v i 2 (x,y) 



where 7 « = 73 = ff±| and ^ = 7l = 
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Proof of Proposition \9.b\ We begin with v 3 since it is more straightforward. Ob- 
serve that for any 7 > 1, using equation (|13|) and the positivity of v 3 , one has 



(Lv 3 )(x,y) = (Av 3 )(x,y) - y 2 d x v 3 (x,y) + a x d xx v 3 (x,y) 

cix s+1 - y 2 \d x v 3 (x, y) \ - cr x \d xx v 3 (x, y)\ 



< -v1{x,y) 



< — mvj (x, y) 



where we define 



m = inf 



C\X 



<5+l 



v 3 {x,y) 



y 2 \d x v 3 (x,y)\ - a x \ d xx v 3 (x, y)\ 



v 3 (x,y) 



If for some choice of 7 > 1 and R > 0, one has m > 0, then it is clear that v 3 
is a local super Lyapunov function. To prove that such 7 and R exist, we use 
the scaling and continuity properties of v 3 which were proven in Proposition 19.31 
Observe that every point {x,y) G lZ 3 (a) can be mapped back to a point (2a, b), 
where (x, y) = S^(2a, b), i — j-, and b = \fty e [—1, 1]. Therefore v 3 satisfies the 
following scaling relations: 



v 3 (x,y) = e s v 3 {2a,b) 



(d x v 3 )(x,y) ^ i 6 '\d x v 3 )(2a,b) 
(d xx v 3 )(x,y) =£ S - 2 (d xx v 3 )(2a,b). 



These scaling relations lead us to choose 7 = ls+t ' w hich is the ratio of the expo- 
nents of £ in Av 3 and v 3 . With this choice of 7, we obtain 



C\X 



.5+1 



y 2 \d x v 3 (x, y)\ - <r x \d xx v 3 (x,y)\ 



v 3 (x,y) 



ci 



(2a) S+1 -e~ 3 (b 2 \d x v 3 (2 ai b)\ +a x \d xx v 3 (2a,b)\) 



v2(2a,b) 



Hence if we define £ 
M 



inf{a;/(2a) : (x,y) € K 3 (a) Ci B C R (0)} and 



sup b 2 1 d x v 3 (2a, b)\ + a x \d xx v 3 {2a, b)\ +v 3 (2a,b). 
fee[-i,i] 

the preceding estimate and the strict positivity of v 3 imply that 

Cl {2a) S+1 -£~ 3 M 



m > 



Since v 3 is C 2 , we know that M < 00 (the supremum is over a compact set). Fur- 
thermore, observe that £* — > 00 as R — > 00. Combining these last two observations 
with the above estimate, we can choose R sufficiently large to ensure that 

Cl (2a) S+1 -£- 3 M > 

and hence that m > 0. We define R 3 and 73 to be the above choices of R and 7, 
respectively, which are valid in lZ 3 (a). Substituting these values in the expression 
for m, we obtain m 3 . This completes the proof of the local super Lyapunov property 
for v 3 . 
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We now turn to proving the corresponding property for v 2 . We start as we did 
for U3, by noting that for any 7 > 1 



(Lv 2 )(x, y) = -h(x, y) + (a x d xx v 2 + a- y d yy v 2 ) (x, y) 
~h - a x \d xx v 2 \ - <Jy\dy V v 2 \~ 



< 



-v 2 (x,y) 



(x,y) < -mv 2 1 {x,y) 



h - <T X \d XX V 2 \ - CTy\dyyV 2 \ 



where in this case we define 
m = inf 

\{x,y)\>R 

As before, we need to show that there exist 7 > 1 and R > such that m > 0. Since 
lZ 2 (a) has two different natural scalings, we decompose this region and handle each 



subrcgion separately Recall the definition of IZ^ and 1Z 



,(2) 



that lZ^\a) scales well under S^' and lZ y 2 '(a) under S 1 



(i) 



,(2), 



: (2) 



,W - 



inf 

\(x,y)\>R 



h - <r x \d xx v 2 \ - <Ty\d yy v 2 \ 



from (|23p and observe 
We define as 



(2) 

We begin with 7Z 2 (ce) since the analysis in this subregion is very similar to the 
previous analysis for V3. 

First, note that the circle of radius r = 2(a 2 + 1) centered at the origin is 

(2) 

contained in A(a, 1, 1). Hence any point in 72-2 (at) can be connected by a radial 

(2) (2} 
line contained in 7Z 2 (a) to the part of this circle contained in 7Z 2 (a). It follows 

(2\ (2) 

from this that any (x, y) € 7Z 2 (a) can be written in the form (x,y) = (a,b) 
where £ — and (a, b) is a point on the circle of radius r centered at the 

origin. Therefore, 

V2(x,y)=£ s v 2 (a,b) (d xx v 2 )(x, y) - t 5 ~ 2 (d xx v 2 )(a, b) 

h(x, y) = e s+1 h(a, b) (d yy v 2 )(x, y) = t s ~ 2 '(d yy v 2 )(a, b) . 

Again, by analogy to the previous case, these scaling relations lead us to choose 

7 = , which is the ratio of the exponents of I in Tv 2 and v 2 . With this choice 

(2) 

of 7, if (x, y) = S\ (a, b), we have that 



h - a x \d xx v 2 \ - <7y\dyyV 2 \ 



(x,y) 



h — i 3 cr x \d xx v 2 \ - £ 3 a y \d yy v 2 \ 



(a, b) 



Setting 71 



(2) 



{(a,b) G Tl 2 2 \a) :\(a,b)\ = r}, we define 



h(a,b) 
(a,b)G7l< 2) v1(a,b) 



inf 



and M 



sup 

(a,b)ea ( 2 2) 



(T x \d xx V 2 \ + <7 y \d yy V 2 \ 



v. 



7 



(a, b), 



Letting = M{\(x,y)\/r : (x,y) G TZ 2 2) (a) n B C R (0)} = R/r, we observe that 
> p — l~ 3 M. Because h and v 2 are C 2 in (x,y) and strictly positive and 
TZ 2 is a compact set, we conclude that p > and M < 00. Hence one can choose 
R sufficiently large in order to ensure that > p — £^ 3 M > 0. Again, we 



denote these specific choices of R and 7, which are valid in lZ 2 2 \a), by R 2 Z> and 



?(2) 



72 



(2) 



Substituting these values in the expression for m( 2 \ we obtain 



,(2) 



This 



PROPAGATING LYAPUNOV FUNCTIONS 



21 



completes the proof that v 2 is a local super Lyapunov property in the subregion 

We now turn to region TZ 2 X \a). Every point (x,y) € H^\a) can be mapped 
back to a point (a, b) on the curve {ab = a} such that (x,y) = S^ 1 \a,b), where 

2 

£ = ^-^j 3 , a = a^(xy 2 )^, and b — a~^(xy 2 )^ . Hence we get the scaling relations 

v 2 (x iy ,l) = £ S v 2 (a,b,£- 3 ) (d xx v 2 )(x iy ,l) = £ S - 2 (d xx v 2 )(a,b,r 3 ) 

h{x,y,l)=£ S+1 h(a,b,r 3 ) {d yy v 2 ){x,y,l) = £ S+1 (d yv v 2 )(a,b,£- 3 ) . 
Now using the scaling map S^p to map (a, b) ^ (a, 1) we obtain 

v 2 (a,b,£- 3 ) = b 5 v 2 (a,l,r 3 ) {d xx v 2 )(a, b, £~ 3 ) = b s - 2 {d xx v 2 ){a, 1, £~ 3 ) 
h(a, b, r 3 ) = b s+1 h(a, l,r 3 ) (d yy v 2 )(a, b, £~ 3 ) = b 5 - 2 (d yy v 2 )(a, 1,£~ 3 ) . 

Again, these scaling relations in 1Z 2 (a) lead us to choose 7 = ||^| , which is the 
ratio of the exponents of £ in Tv 2 and v 2 . Combining these two sets of scaling 
estimates and setting 7 = 5(1 — 7) + 1, we get that 

-h - <7 X \8 XX V 2 \ — CFy\dyyV 2 \ 



= w 



{x,y, l) 

h - (b£)~ 3 a x \d xx v 2 \ - b~ 3 a y \d yy v 2 \ 



(a,i,r 3 ) 



±(l-(M)- 3 * x ^-b- 3 * v ^)](a,l,l 



dyyV 2 \ 



We have organized this calculation a bit differently for reasons which will become 
clear momentarily. Based on the preceding calculation, we define 



/ \ \ . t h(a,l,X') 
p{\) = inf -j- — —— 
A'e[o,A] vj(a, 1, A') 



and Mi (A) = inf 



\d xx v 2 (a, 1, A')| 



A'e[o,A] h(a,l,\') 



and M 2 (A) = inf 



\d yy v 2 {a,l,X')\ 



A'e[o,A] h(a,l,\') 

Notice that, unlike the previous calculations, we have made the constants p, Mi, 
and M 2 depend on the maximal value of A. This is because in our current setting 
we require more precise information about these constants than merely that they 
are finite and positive. 

We set £* = inf{(x/ay) 3 : (x,y) € 1Z 2 (a) n B R (0)} and we observe that since 
6>1, 

m (1) > p(C 3 )(l - CV x Mi(C 3 ) - (JyM 2 (£- 3 )) . 

We wish to conclude that the right-hand side of the above expression is strictly 
positive. To conclude this, however, we need to understand the behavior of Mi (A) 
and M 2 (A) as A — > 0. By direct calculation from the explicit formula for v 2 , we see 
that 



Mi(0) 



2(5(2(5+1) 



and M 2 (0) 



(<S +!)(<* + 2) 



a° a 
Since Mi (A) and M 2 {\) are both continuous functions of A on (0, 1] with finite limits 
as A — > 0, and since 1* can be made arbitrarily large by choosing R sufficiently large, 
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for any e > we can choose R large enough to ensure 

1 - C 3 ^M l(C 3 ) - ^M 2 (C 3 ) > 1 - ^ S + f S + V - e . 

We conclude that as long as a y( s+1 J( s + 2 ') < \ j we can always choose i? large enough 
to guarantee that is positive. This last inequality holds whenever a is suf- 
ficiently large. Again, we denote these specific choices of R and 7 , which are 
valid in 1Z 2 1 \a), by R^ and j^K Substituting these values in the expression for 
toW, we obtain m^. Choosing m 2 = minjm^ 1 ', m^}, 72 = min^ 1 ^, J^}, and 
i? 2 = maxji?^ 1 ', -R^} completes the proof that v 2 is a local super Lyapunov func- 
tion in the entire region lZ 2 (a). The more detailed statements of the behavior in 
and merely serve to summarize the above points. □ 

10. Construction of a Global Super Lyapunov Function 

We now patch together the three local Lyapunov functions that are defined in 
distinct regions of the plane in order to produce one smooth, global Lyapunov 
function defined on the entire plane. To do this, we use the standard mollifier <j)(t), 
a smooth, increasing function which varies from to 1 and is suitably normalized to 
integrate to unity on the whole real line. Specifically, we take 4>(t) = ^ j_ x ip(s)ds 
with m = tp(s)ds where 




for < t < 1 
otherwise . 

Next, we define the functions hi(x 7 y) and h 2 (x,y) as follows: 

x a 
The function h\{x,y) — on one boundary of the wedge-shaped region TZ\{a) n 
7?.2(a); hi(x,y) = 1 on the other boundary of this region; and hi varies smoothly 
between and 1 in the interior. Similarly, h 2 {x,y) — on one boundary of the 
funnel- like region 72-2 (a) H 72.3(a) and h 2 (x,y) = 1 on the other boundary. Thus, 
outside of a fixed ball, we define our global Lyapunov function V to agree with 
the local Lyapunov functions in subregions of their domains of definition and to be 
a smooth, convex combination of the two local Lyapunov functions in regions of 
intersection. In particular, let V(x, y) be given by 

'v 1 (x,y) for (x,y) G 1li(a) n TZ 2 (a) c 

Vi(x,y) for (x, y) G 72 1 (a) H IZ2 (a) 

v 2 {x,y) for (x,y) G K 2 (a) n 72i(a) c n TZ 3 (a) c 

V 2 (x,y) for (x,y) G TZ 2 (a) DTZ 3 (a) 

v 3 (x,y) for (x, y) G 1Z 3 {a) n 7Z 2 (a) c 

otherwise 

where 

Vi{x,y) = [1 - (j){h l {x,y))]v 2 {x,y) + ^{h^x.y^v^x^y) 
V 2 (x,y) = [1 - (f>(h 2 (x,y))]v 2 (x,y) + 4>(h 2 {x,y))v 3 (x,y) . 



V{x,y) = < 
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We then choose V(x, y) € C 2 (R 2 ) to satisfy 

jv(x,y) iorx 2 + y 2 >p 2 

V[x, y) = < , 909 

I arbitrary positive and smooth for ar + y < p 

where p will be specified below. 



Remark 10.1. At the start of the Lyapunov construction in Section \ 9.1l we fix 
a choice of S > when defining v\ . This choice is then propagated through our 
construction and is explicitly present in the definition of vi and V3 . During the 
analysis ofv^, we note in Provosition \9.S\ and Provosition \9.3\ that we must choose 
S 6 (0,f)- Except for this one restriction, we are free to choose 6. Hence our 
construction of V depends on two parameters 6 and p. As we will summarize in 
Proposition \10.2\ below, 5 gives the power of the polynomial growth in all but the 
pure, positive x-direction. On the other hand, p sets the distance from the origin 
past which the Lyapunov estimates are valid. 

Consolidating our results on the scaling behavior of the functions Vi, i = 1, 2, 3, 
we obtain the following: 

Proposition 10.2. Fixing a 8 € (0, |), there exists positive constants c and C , so 
that c\(x,y)\ s < V(x,y) < C\(x,y)\ S = C\(x,y)\M . 

Proof of Provosition \10.2[ After observing that 6 < 5 = |<5+ 1 on (0, |), the result 
follows quickly from Proposition 19.31 and the definition of v\ from Section 19.11 On 
the right half-plane, the result follows from the definition of v-y. On the left half- 
plane, one can either use the the scaling relations or the explicit representations 
given in (|21[) and (|20| to obtain the desired inequalities. □ 

The following proposition states that V is a super Lyapunov function. Therefore, 
one of our main theorems, Theorem [T] from Section[3j is an immediate consequence 
of this proposition. 

Proposition 10.3. For any 6 G (0, |), there exists a p from the definition ofV so 
that V(x,y) is a global super Lyapunov function on M 2 . 

In light of Proposition 19.11 and Proposition 19. 61 the main missing ingredient in 
the proof of Proposition 1 1 . 3l is the verification that V is a local Lyapunov function 
in the patching regions. This is the content of the next proposition; we prove it 
before returning to the proof of Proposition 110.31 at the end of the section. 

Lemma 10.4. For any S £ (0, |), Vi(x,y) is a local super Lyapunov function on 
lZi(a) H lZ2(c < ) and ^{x, y) is a local super Lyapunov function on lZ2{ a ) H TZs(a). 

Proof of Lemma \10.4\ Let m„ Ri and ji for i — 1,2,3 be the constants from 
Proposition 19.11 and Proposition 19.61 Next define to* = min{mi, TO2, to-3}, R* — 
max{i?i, i? 2 , R3} and 7* = min{7i, 72, 73} — ff^f ■ We further increase i?* so that 
if (x,y) £ TZi and |(x, y)| > i?* then Vi(x,y) > 1. In proving that V2 is a local 
super Lyapunov function we need to recall the more refined version of the super 
Lyapunov property in IZ2 given in Proposition 19.61 

We address each of the claims in the lemma separately. We begin with the proof 
that V\ is super Lyapunov since it is the most straightforward. If p > R* , we have 



24 



AVANTI ATHREYA, TIFFANY KOLBA, AND JONATHAN C. MATTINGLY 



that for all (x,y) € fti(a) n TZ 2 (a) n B°(0) 

(LVi)(a;,y) = (1 - </>(hi(x, y))){Lv 2 ){x, y) + <f>(hi(x, y))(Lvi)(x, y) + E^x, y) 

< -m»[(l - ^{h^x^y^vfix^y) + 4>{hi{x, y^vj 1 (x, y)] +E 1 {x,y) 

< -m*[Vi(a;,t/)] 71 + Ei(x,y) by convexity 

< -m,(l-Mi)[^i(a;,y)P 
where Mi and y) are defined as 

Ei(x,y) 



Mi = sup 

(i,t)£Ki(o)nKi(a) 
\{x,y)\>p 



m*[Vi{x,y)] 



71 



and 



Ei(x,y) = L[(j>{hi{x,y))]{vi{x,y) - u 2 (z,y)) 

+ 2^ — [</>(/n(a;,y))] — [ui(a;,y) - w 2 (a;,y)] 

+ 2 °y [<P(hi {x, y))} ^ [vi (x, y) - v 2 (x, y)] . 

If we can choose p sufficiently large so that M\ < 1, then V\{x,y) will be a local 
super Lyapunov function on 1Z\(a) n 1Z 2 (a). To show Mi < 1, we use the scaling 
properties of v\ and v 2 to map back to a circular arc of radius r = Vet 2 + 4. Let 



\A 2 + y 2 



a = t, 



and b — — . 



(x,y)=sf\a,b) with l>\ 



Then 

(x,y)e^i(a)n^ 2 (a)n^(0) = 

Note that h\{x,y) = h\(a,b), so 

V 1 (x,y) = 1 & V X (a, b) and \Vi (x, y)p = £ s+1 [Vi (a, 6)p . 

As a consequence of these scaling relations, we get that for all (x,y) € TZ\(a) n 
Tl 2 (a)nB c p (0), 

Ek(x, y) = e^afiih^a, b))\b\ (l + ^) [«i(o, b) - w 2 (a, b)] 



+ e s - 1 a(f>"(h 1 {a,b)) 



-a x \b\ a y sgn(b)\ 

+ ) [ui [a, b) - v 2 (a,b)\ 



+ Wfai (a, b))^P- [m (a, b) - W2(a, 6)] 

+ ^ W(M«, b ))^^^-[ Vl (a, b) - v 2 (a, b)]) 

+ l^^{h x (a, 6)) 2fW(6) £ [ Vl (a, 6) - v 2 (a, 6)] . 
a Cy 



Hence we have that 

(24) Mi < sup 

(a,b)eni(a)nK 2 (a) 
\(x,y)\<r 



ei,i(a, b) ei ;2 (a, b) 

m.[Vi(a,&)]Ti £ 2 m4yi(a,6)] 
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where 

ei,i(a, b) = a<p'(hx(a, b))\b\ (l + - 2 ) [«i(o, b) - v 2 (a, b)} 

e h2 (a, b) = acj>"{h x {a, 6))(=^ + ^1^1) [ Vl (a, 6) - *, 2 (a, 6)] 

V ar a J 

+ acp'ih^a, b))^^{ Vl (a, b) - v 2 (a, b)} 



+ a^(/n(o, &))_?j!M JL^fo 6) ~ «a(o, 6)] 



+ a^(/n (o, 6)) 2fTySgn(&) £ [«! (a, 6) - v 2 (a, b)} . 
a dy 

By explicit computation with v\ and v 2 , we can show that ei j i(a, b) is always 
negative for (a, 6) in the desired region. The second term of the sum in (l/T^tl) . 
the upper bound for Mi, can then be made arbitrarily small by choosing £ large 
enough; this corresponds to choosing p sufficiently large. This establishes that 
Mi < f , which completes the proof of the lemma. 

We now turn to proving that V 2 is super Lyapunov. If p > R* , then for all 

(x,y)en 2 (a)nn 3 (a)nB c p (o) 

(LV 2 ){x,y) = (1 - <j)(h 2 (x,y)))(Lv 2 )(x,y) + ^(h 2 (x,y))(Lv 3 )(x,y) + E 2 (x,y) 

< -m»[(l - 4>{h 2 (x,y)))v 2 ,3 (x,y)+(t){h 2 {x,y))vf(x,y)] + E 2 (x,y) 

< — [V 2 (x, y)] 73 + E 2 (x, y) by convexity 

< -m*(l-M 2 )[V 2 (x,y)}^ 



where 



and 



M 2 = sup 



(i,9)eRs(a)nR 3 (a)nBj(0) L TO * [^2 (x, y)p 3 



E 2 {x,y) 



E 2 (x,y) = L[0(/i 2 (£,y))](u3(>,y) - U2(af,!/)) 

+ 2vx-^[<j>(h 2 (x,y))] — [v 3 {x,y) - v 2 {x,y)\ 

d d 
+ ^y-^[<t>{h 2 (x,y))] — [v 3 (x,y) - v 2 (x,y)] . 

If M 2 < 1, then V^x, y) will be a super Lyapunov function on 1Z 2 (a) n 7^3 (a). To 
show M 2 < 1, we use the scaling properties of v 2 and v 3 to map back to a vertical 
line. Let 

£=—, a = 2a, and b = yVt ■ 
2a 

Then 

(x,y)£ll 2 (a)nll 3 (a)nB c p (0) =>■ (a:,!/,l) = f (a,M" 3 ) 

where |6| £ [^=, l] and £> 1. Note that ^(x, y) = /i2(a, so V 2 (x, y) = V 2 (x, y, 1) 
satisfies 

V 2 {x,y,l)=t*V 2 (a,b,r 3 ) and {V 2 (x,y,l)}^ = £ S+1 {V 2 (a,b,r 3 )}^ 

where 

V2O, y, A) = [1 — <j){h 2 (x, y))]v 2 {x, y, A) + <j){h 2 {x, y))v 3 (x, y) . 
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Now we have that for all (x, y) £ TZ 2 {a) (~l TZ 3 (a) n B c (0), 



E 2 {x,y) = t 



6 + 



i</>'(h 2 (a,b)) 



a 



(-5a V - 2aa y )[v 3 {a,b) - v 2 (a, b,r 3 )} 







e s + i 0'(h2(a,b)) „ _ V2iaAe - 3)] 

a ay 

£ S +i r(h2{a > b)) (-2aba y )[v 3 (a, b) - v 2 (a, b,£~ 3 )] 



+ t 



+ £ 



a 

s^ 2 <j>'(h 2 (a,b)) 



(-b 2 a x )[v 3 (a, b) -v 2 {aA^)] 



6 4 h(a,&)- W2 (a,6,r J )] 



c $_ 2 <t>' {h 2 (a, b)) 2 d 3 
I (-2b a x ) — [v 3 (a,b)-v 2 (a,b,£ )\ . 



Dehne N(X*) as follows: 
(25) N(X*) 



sup 

Ae(o,A*] 



e 2 ,i(a, b, A) 



e 2 , 2 (a, 6, A) 



m,[V 2 {a,b,X)}^ sTlm*[V 2 {a,b,\)Yi* 



where 



e 2 ,i(a, b, A) 



e 2j2 (a, b, A) 



^ {h2{a ' b) \ -5a 2 b 2 - 2aa y )[v 3 (a, b) - v 2 (a, b, A)] 

<t>\Ma,b)) f _ inhrT) d_ 
dy 



-(~^aba y )^[v 3 (a,b) - v 2 (a,b, A)] , 



0"(h 2 (a,b)) 



(-2aba y )[v 3 (a, b) - v 2 (a, b, A)] 



<//'(h 2 {a,b)) [ 
<f>'(h 2 (a,b)) u4 



-b 2 a x )[v 3 (a, b) - v 2 (a, b, A)] 



b [v 3 (a,b) - v 2 (a,b,X)} 

+ <t>\h 2 {a,b)) d 
a ox 

Note that for any A* > 0, we can choose p sufficiently large to force M 2 (which, 
we recall, depends on p) to be less than N{\*). Ultimately, we will choose A* 
sufficiently small so that iV(A*) < 1. The second term of the sum in (1251) can be 
made arbitrarily small by increasing the size of £; again, increasing the size of I 
corresponds to increasing p. We now address the first term of the sum in (|25[) . 
From Lemma flO.51 which is stated and proven bellow, we see that we can chose the 
parameters to make this term negative. 

Combining all of these observations, we have demonstrated that M 2 < 1, which 
completes the proof of the lemma. □ 



Lemma 10.5. There exist positive constants c\ and c 2 in the definition of the 
Poisson equation for v 3 (x,y), and positive a and A* such that for all A £ [0, A*], 
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the following inequalities hold for a = 2a and \b\ 6 ["T^ji l] •' 
(26) t&(a,&) -v 2 (a,b,X) > 



(27) 



^M)-^(a,M) 



> 



Proof of Lemma \10.5[ Recall that, from ([20]) . vs(x,y) can be represented as 



(28) 



f3(z,y) = x 



{f + c 2 )E (x , y) | 



„<5rl 



where t = inf{f > : \Zt\ £ [— \/2a, V2a]} and Z t is the process given in (IT51) . 
Note that the expectation in ([28)1 can be written as the solution to a second-order 
ODE, namely: 

E (x^)[e 5r ] = 9e{Vex y) 
where g t {z) solves the following boundary value problem with e 



2a • 



(29) 



(ecr y g'J(z) + \zg' e (z) + Sg e {z) = for z G (-1, 1) 
\g e (-l)=g e (l) = l. 



Define go(z) to be the solution to the limiting ODE in (f29| when e = and note 
that go(z) can be computed exactly for initial conditions z =/= 0: 

(30) .g (z) = — j-r. 
Now, let vo(x, y) be defined as 

(31) V {X: y )=X S\^ + c ^g o{V ^ y) _^ 



We address the first difference between V3 and v 2 in (126)1 as follows. We write 

(32) v 3 (a, b) - v 2 (a, b) =u 3 (a, b) - v° (a, b) 

(33) + v° 3 {a,b)-v 2 {a,b,0) 

(34) +W2(a,6,0)-«3(a,6,A). 

To show that this difference is positive, we will show that v®(a, b) — 1*2(0, b, 0) > 
and that the other two differences are small in comparison. Similarly, for the 
difference between the y-derivatives of V3 and v 2 in (|27[) . we write 



(35) 
(36) 
(37) 



dv 3 , . 9u2. . 
^(a,6)-— (a,6,A) 



=6 



— (a,&)-— (a,6) 



+ 6 
+ 6 



and again, we will show that &[^-(a, 6) — ^-(a, 6,0)] > and the other two 
differences are small in comparison. Specifically, we demonstrate that there exist 
positive constants c\ and c 2 in the Poisson equation for 1*3 such that the differences 
in (f3"3"|) and (|3l)|) are positive; and then, that there exists a sufficiently large such 
that the differences on the righthand sides of (|3"2"j) and ([3"5"j) are comparatively small; 
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and last, that there exists a A* such that (|34|) and (|37|) are comparatively small for 
all A € [0, A*]. For the first of these claims, note that 

(2af^ 



(38) 



v°(a,b) - v 2 (a,b,0) = 



_( a ,6)-— (a,6,0) 



|6| 



\b\ 



5+ 



where q and q are given by 



g(6) = 2*'+* (f + c 2 )H« 



Cl 



l&l 



.5+1 



q(b) = -(6 + |)2* 5 +i (f + c 2 ) |6|f + 5 + 1 



and ci 



-r and C2 = 



with a so that v 2 and 1/3 have identical scaling in a. Moreover, as we demonstrate 
below, c\ and c~ 2 can be chosen independently of a. It is clear that q is a monotone 
decreasing function of \b\; hence it is minimized at the right endpoint of the interval 
for |&|, that is, |b| = 1. Thus if we can show q(l) > 0, then it follows that 



We note that c\ and c 2 are chosen to scale 



(39) 



*i(a,6)-^M,0) 
dy dy 



> 



for all |6| € [2-^,1]. If we can also show q(2 2) > 0, then from (j3"8")l . we conclude 
that 

u§(a,2-5) _ t; 2 ( a ,2~5,0) > 0. 

Combining this with (j59")l gives the desired positivity of (j3"3")l on the whole interval 
\b\ e [2-2,1]. Hence, we need only verify that there exist positive values of 6\ and 
C2 such that 



(40) 
(41) 



(<| +c -)2-*_f2-^+5) 



- 1 > 



g(l) = -(<5 + f )2* a +* (f + c 2 ) + 6 + 1 > . 

The verification of this is elementary and we omit the details. 

We remark that c\ and c~2 can be chosen independently of a, since the above 
inequalities have no dependence on a. Thus, choosing positive 5\ and c~2 such that 
(|40| and flJT]) are both satisfied, we get that for all \b\ <E [2 _ 5, 1], 

v° 3 (a,b) -v 2 {a,b,0) > a 25+1 2i 5+ iq{2-^) > 
dv 2 



dy 



(a, 6,0) >a 25+1 2 25+1 q(l) >0. 



We turn our attention to making the differences 

1/3(0, 6) — 1/3(0, 6) and 6 

comparatively small. Note that 

,(a, 6) - v° 3 (a, 6) = a 25+1 2§ A "+i (f + c 2 ) [g e (b) - 5o (&)] 



dv 3 dvl 
— (a, 6) -—(a, 6) 



«3l 



— (a, 6) -—(a, 6) 



a 25+i 2 § 5 +f6(f 



c 2 )[^(&)-5oW] 
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To be precise, we will show that 



\v 3 (a,b)-v° 3 (a,b)\ < § 



Z(2"*)] 



dv 3 dvl 
— (a, 6) -_ (a> 6) 

This is equivalent to establishing that 



<*[« 



3(1)] 



(42) 
(43) 



„2<5H 



1 2S*+S| 



c 2 )[g e {b)-9o{b)] <\[a 2S+ ^+i q 



a 



2*+l f 



< | | a 



(2-i)] 



Observe that the same powers of a appear on both sides of each of the above 
inequalities. Therefore, to prove (|42[) and (|43|) - it suffices to show that g e (b) and 
(7^(6) converge uniformly to go(b) and g' (b), respectively, for \b\ G [2 _ 5 j l] as e = 

tt > 0. Both of these uniform convergences follow from classical results; see, for 

example, [AS92J . 

Since e = ^- , we can choose a sufficiently large to guarantee that both (|4"2"j) and 
(l43l) hold and that v-i remains a local super Lyapunov function on 7^.2 (en) (recall 
that in Proposition 19.61 a lower bound on the size of a was imposed) . Finally, by 
choosing A* sufficiently small, the differences in (f34| and (|37| can be made small 
for all A G [0, A*]. This is an immediate consequence of the fact that V2(a, 6, A) is a 
C 2 function of A G [0,1]. □ 

Having established the super Lyapunov property in the patching region, we re- 
turn to the proof of the main result of this section. 

Proof of Proposition \10.3[ The local super Lyapunov condition has now been proven 
in regions; namely, for v\ in Proposition 19. 11 for V2 and v 3 in Proposition 19. 6[ and 
for the patched functions V\ and V2 in Proposition ll0.4l All that remains is to make 
a global choice of constants. The constant p from Proposition 110.41 was chosen to 
be valid in all regions. It is sufficient to choose 

M = min {m st (l — Mi),m*(l — M 2 )} < to*, 

b = sup{\{LV){x,y)\ :x 2 +y 2 <p 2 }, and 

, 5(5 + 5 
7 = mm{7i , 72 , 73} = — - - . 

These choices guarantee that for all (x, y) G R 2 

(LV)(x,y)<-M[V(x,y)r + b. 

□ 



11. Existence and Positivity of Transition Density 

Having established the existence of a global super Lyapunov function, we now 
make a small detour to prove the existence of a smooth density with appropriate 
positivity properties. These results provide the missing ingredient to prove the 
ergodic result stated in Theorem [5J namely a minorization condition. It is worth 
noting, however, that proving the minorization condition is not our sole goal. In- 
deed, if it were, we would not need all of the results of this section: we could simply 
use smoothness and appropriate open set controllability results. See for example 
[MSH02IIMMP09] . 
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Instead, our interest is larger, motivated by two concerns. First, we wish to 
understand the general structure of the invariant measure, not merely its unique- 
ness. Second, we want to take this simple example to highlight some techniques, 
different than those often used, which can be applied in more general situations to 
address questions of positivity. We feel that these methods convey more intuition 
and better allow for the inclusion of a priori facts about the dynamics. 

11.1. Positivity when a x > 0. When a x > (since we always assume cr y > 0), 
the system is uniformly elliptic, and everything is relatively straightforward. Since 
the diffusion associated with ((T|) has a constant, positive definite diffusion matrix, 
classical results guarantee the existence of a function p : (0, oo) x R 2 x R 2 — > (0, oo) 
such that p is jointly continuous, p t {z,z') is strictly positive for all (t, z,z'), and 
such that for all measurable sets A 



We summarize these results for future reference in the following Proposition. 

Proposition 11.1. If a y > and o~ x > then for all t > 0, P t has an everywhere 
positive density p t (z, z') with respect to Lebesgue measure which is smooth in both 
z and z' . 

11.2. Positivity when a x = 0. When a x = (but a y still positive), the situa- 
tion is more delicate. We begin by observing that the generator of the associated 
diffusion is still hypoelliptic; to see this, we write the generator of the diffusion as 



Observe that [[X, d y ),d y ] — —2 and hence the relevant ideal in the algebra generated 
by X and d y is of full rank. In turn, this ensures the existence of a continuous 
function p so that (j44j) holds. The main difference between this and the setting of 
Section ril.il is that it is no longer immediate that pt(z, z') is positive for all t > 
and z, z' e K 2 . In fact it is not true. 

Intuitively, it is clear that if z is in the left-half plane and z' in the right-half 
plane then pt(z, z') should be zero, since there is no way to move across the y-axis. 
On the other hand, it is reasonable to expect that given any z' in the left-half plane, 
there exists a T = T(z') such that p t (z, z') > for alH > T and z £ K 2 . 

There are a number of ways to prove such a result. The most generally applicable 
and powerful technique is to leverage geometric control theory to show that the 
support of Pt(z, ■ ) contains a sufficiently large, bounded region of the left-half 
plane for any z and t sufficiently large. From this, for example, one can show that 
(Xt,Yt) is sufficiently smooth in the Malliavin sense (which it is), and deduce that 
Pt(z, z') is strictly positive in the interior of the support. 

Alternatively, one can use sufficiently quantitative open-set controllability results 
to extend the very local positivity which follows just from the joint-continuity of 
(z,z') i-> Pt{z,z'). This is the method employed in [MSH02, MMP09 in various 



Here we take an approach most consonant with the first option. However, rather 
than merely citing the appropriate geometric control theory results, we construct 
an explicit series of simple controls to prove the positivity condition we require. 



(44) 




L = X + -a y dl . 



forms. 
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The subsequent discussion is lengthier, but we feel that it is more intuitive and is a 
useful complement to more general control theoretic arguments, especially for the 
uninitiated. 

Before turning to the existence of a positive density for (fTJ) , we first consider an 
analogous calculation in a simpler setting. Consider a smooth map cj>: M. m —> M. d 
where m > d and let T be a non-degenerate Gaussian probability measure on W n . 
Consider the push forward of T by 4>, denoted by Ftp' 1 and defined by Tcj)~ 1 (A) = 
T(<j)^ 1 (A)). For the measure Ycf)^ 1 to be absolutely continuous with respect to 
Lebesgue measure, it is necessary and sufficient that 



(See [Bel95, B el04| for more details.) Supposing that this condition holds, we let 7 
denote the density of with respect to Lebesgue measure. We are interested in 
when 7(2) is positive at a given point z E M. d . It is a simple exercise in calculus to 
see that j(z) > if and only if there exists a x G R m with 7(2;) = z and for which 
(D<t>(x))(D(f) T (x)) is a non-degenerate matrix. The first condition ensures that 
there is a way to reach z; that is to say, z is the image of R d under <f>. The second 
ensures that an infinitesimal piece of volume, and hence probability, is brought with 
x when it is mapped by (f>. 

This intuitive fact has a counterpart in stochastic analysis. While these ideas rest 
on the foundation of Malliavin calculus, the closest analogue is found in the work 
of Ben Arous and Leandre [BAL91b , BAL91a and the subsequent presentation by 
Nualart [BN98j . We begin by identifying the map in question. 



To any U e L 2 ([0, T], R d ), we associate {(A"f , Y t u ) : t € [0,T]} which solves the 
system of equations 



In the following discussion, we will refer to U as the control and denote by (Xf, Y d ) 
the solution to the deterministic system of equations (|4|); that is, the system ([45]) 
with U(x, y) identically zero. It is also worth mentioning that for U € L 2 ([0, T],M), 
t i-» (Xf , Y t u ) is continuous on [0, T}. 

In analogy to the discussion at the start of the section, for T > and z € K 2 , 
we will consider the map $ T , Z : L 2 ((0,T];IR) R 2 defined by U ^ ($W ,$^ 2) J = 
(Xj,, Yjl) and (Xq , Fj 7 ) = z. Translating [BN98] to our current setting, we obtain 
the following theorem: 

Theorem 4. px(z,z') > if and only if there exists a U £ L 2 ([0,T],M) so that 
&t,z(U) — z' and furthermore the matrix 



is nondegenerate. Here D represents the Frechet derivative. 

The condition from |B AL91bl IB AL91al IBN98] that the SDE under consideration 
have coefficients which are bounded with all derivatives bounded can be removed 
by a standard localization argument. The key step is knowing that |(A"t,Yt)| is 
almost surely finite which follows from the Lyapunov function we constructed in 
the preceding sections. 



T{x 6 M m : Det\(D(i>)(x){D(j) T )(x)\ = 0} = . 



(45) 




(Af ) 2 - {Y t u f 
2X^Y t u + U t 



(46) 
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This matrix M is just the product of the Jacobians introduced in the motivating 
discussion at the start of this section, translated to our current setting. To facilitate 
calculations and to better see this connection, it is useful to observe that for any 

t) e K 2 

(47) (M T , z (U)r,,r ] )= [ {J^ T e 2lV ) 2 ds 

Jo 

where e 2 = (0, 1) and J^ t is the Jacobi flow (the linearization of the SDE/controlled 
ODE) defined by 



J. 



u 



( d^l(U) d<S>[]l(U) \ 
d^iiU) d<S>i 2 l(U) 



\8*$(17) d<S> {2 l(U) / 

We now state a simple condition which ensures the nondegeneracy of M. It 
captures the fact that as long as there is some twist in J SjS + e then J SiS+c e 2 and e 2 
will not be co-linear, and hence ( J SjS + e e2, rj) + (e 2 , rj) ^ for any r\ ^ 0. Combining 
this with the continuity of s H> J s ,t produces the desired nondegeneracy of (|47| . 
The following lemma follows this outline, providing a condition which ensure such 
that the system has the desired twist. 

Proposition 11.2. To ensure the nondegeneracy of Mt, z {U), it is sufficient that 
there exist a t £ [0, T] so that &t ,z(U) =/= 0. 

Proof of Provosition \11.2L Since 1 i-> $ tiZ (U) is continuous, we can without loss of 
generality assume that to £ (0,T) and pick an e > so [to — e 7 t ] C (0,T) and 

$t,z( U ) > for a11 * € [to - e,t ]< 

The nondegeneracy of Mt,z(U) is equivalent to 

inf (Mt z(U)r),r)) > 0. 



In light of (|47|) . we see that 

(M T , z (U)r,,r,) > [° (J^ T e 2 ,r,fdt = f° ( J t ^ e 2 , (J^ 0iT )*v) 2 dt 

where (J^ g T )* is the adjoint of the matrix J^ g T . Combining these two observations, 
we see that for some positive constant c, depending on U, 

rto+e 

inf (M T z (U)n,ri) > c inf / (J^ e 2 ,ri) 2 dt . 

Since (t, rf) i— > (jjf tg e 2 , rj) is jointly continuous, it is sufficient to show that 
(48) ( j£ o e 2) e£) ^ for all t G [to - e, to) 

for some positive e where e^ is the standard choice of vector perpendicular to e 2 . 
Since J^ Q tg e 2 = e 2 , this guarantees that for every given n ^ 0, one has (Jfjo^iW) ^ 
for t in some open interval of [to — e,to]- The continuity in n then ensures the 
infimum over all n with \r)\ — 1 is still positive. 

To establish (j48|) . we appeal directly to the equations. We see that as long as 
®t,z{U) 7^ to t £ [t -e,fo] then {J s ^n,^) ^ for all s,r with to — e < s < r < to 
and all rj ^ 0. This is due to the fact that as long as the trajectory is not at zero, 
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the linearization rotates any vector a nontrivial amount over any interval of time. 
In particular, the linearization satisfies the equation dtJ s ,t — AtJ s ,t for t > s with 
J s . s equal to the identity matrix and 

(AQ\ A-( 2Xt ~ 2Y *\ -91? ( C ° s6t " Sm M 

(49) At ~\2Y t 2X t )- 2Rt [ S me t cosej 

where (X t ,Y t ) = (R t cos t , Rt sin 9 t ). □ 

Remark 11.3. The proof of Proposition Vll.2\ gives very appealing intuition con- 
cerning the positivity of the transition density. Stochastic variation is injected at 
every moment of time in the y-direction. However, the deterministic part of the 
flow is rotating at every point except the origin, as is seen from the calculation in 
(|49[) . This rotation ensures that there is stochastically independent variation in two 
linearly independent directions. 

As a consequence of Proposition lll.2l to invoke Theorem@]to prove the positivity 
of p t {z,z') for two given points z,z' € R 2 , we simply need to find a control U for 
which &t,z(U) = z' and for which the path does not spend all of its time at the 
origin. Since the path is continuous in time, this last condition poses no restriction if 
either z or z' is not the origin. If both z and z' are the origin, it is still elementary to 
find a control satisfying the second condition which still also satisfies &t,z (U) = z' . 

We collect these last observations in the following lemma which combines Propo- 
sition QX2] and Theorem |U 

Corollary 11.4. The transition density p t satisfies pt{%, y) > for a given T > 
and x,y € R 2 if there exist a U G -L 2 ([0,T],R) such that ^>t,x{U) = y and there 
exists an s £ [0,T] so that & s ,x(U) ^ 0. Similarly, pT{x,y) = if there exists no 
U £ L 2 {[0,T],R) such that $t,x(U) = y. 

We now build the required controls. All the controls we design will take the 
form U t = u(X^ ,Y t u \t) for some piecewise smooth u : R 2 x [0,T] — > R. At first 
glance, this might seem an implicit specification which risks being ill-defined, since 
(X^ , Yj 7 ) depends on the function U t we define through ([45|) . However, a moment's 
reflection shows this not to be the case, since in this setting (X^ , Y^ 7 ) can be defined 
in a self-contained way as the solution to the ODE 

Y t u = 2X?Y t u + u(xV,Y t u ,t). 

Then, with this solution in hand, one can define Ut = u(X^ , Yf , t). 

Lemma 11.5. Let z* — with re* < be fixed. There exists a finite > 

such that for all zq — (xq, yo) & R 2 and for all T > T Sf {z^), there exists a control 
U G L 2 ([0,T],R) for which the controlled system (|45|) satisfies 

{Xq,YU) = z , (Xt,Yt) = z* 

and such that Mt, Zo (U) is nondegenerate. 

Proof of Lemma \11.5\ We begin with the case of = (x*,y*) with ^ 0. The 
remaining case will be considered at the end of the proof. The critical feature of 
z* with ^ is that there exists a deterministic orbit which begins on the y-axis 
and flows in finite time to the point z*. 
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{%* , y* ) 




(xo,yo) 



X 



Figure 3. A representative example of a control path used to 
move from any point in R 2 to any point in the left-half plane. 

Let B be the closed ball of radius |a;*|/3 about the origin intersected with the 
negative half-plane {(x, y) : x < 0}. We begin by setting U t = sgn + (y ) for t G [0, 1]. 
(Here sgn + (a;) = x/\x\ for x ^ and 1 if x = 0.) Define T\ to be the first time 
t > 1 such that {Xf,Y^) G B where (X$ ,Yg) = z Q . We will set U t = for 
t G [l,Ti + s] where s > is a parameter we will vary in our construction. By 
driving with \Ut\ = 1 on the time interval [0,1], we have ensured that ^ 0, 
which in turn implies that T\ < T-j* for some finite, zo-independent constant T*. 
(In light of the discussion at the end of Section [3 7\* < 1 + t^t-) Let C* be the 
orbit of the deterministic system which passes through the point but which is not 
contained in the set B. Recall that this orbit is a circle tangent to the origin. We 
now define T 2 to be the infimum over time t > T\ + s such that (X\ V , Y^) 6 C*. If 
we define U t = sgn+ (Y t u )M - 2X^Y t u for t G [T x + s, T 2 ] , then Y t u = sgn+ (Y t u )M, 
and Y t u = Y^ i+s + M sgn+ (Y^ i+S ) (t - T t - s) for t G [T x + s, T 2 ] . Hence for M large 
enough, we can ensure that X t moves very little in the time it takes Y t to grow 
sufficiently to cross C* . This allows us to prove that T 2 — (7\ + s) is bounded from 
above for any sufficiently large fixed choice of M with a bound which is independent 
of s and zq since {X% +s , Yj? +S ) G B. For the same reason, by choosing M large 
enough we can ensure that \X% \ < 2|ac*|/3. Fixing such an M, we define T3 to 
be the infimum of times t > T 2 such that (X^ ,Y t u ) = z*. For t G [T 2 ,T 3 ], we set 
Ut — 0. Since M is fixed and \Xj, \ < 2|x*|/3, clearly T3 — T 2 is bounded uniformly 
for all s > and z G K 2 . If we view T 3 - (7\ + s) as a function of (X^ 1+S ,Y^ 1+S ), 
then it is continuous, since the governing ODEs depend continuously on their initial 
conditions. Since (Xj, +s , Y^ i+S ) G B, which is a compact set, we know that there 
exists a finite bound r* so that T 3 - (T x + s) < r* for all (X% i+a , Y^ i+S ) G B. 
Because {Xj, +s ,Yp, s ) is a continuous function of (zq,s) and T\ an continuous 
function of zq, we see that (zo,s) i-> T 3 is a continuous function. Since T\ and 
T3 — T 2 are bounded uniformly in (zo,s), we conclude that if T : (zo ; s) n- T3, 
then there exists a so that T(zo,0) < for all zq. Since s T(zo,s) is 
continuous and T(zq, s) — ¥ 00 as s — > 00, we conclude that for any t > and 
zq G M 2 there exists a s(t, Zo) so that T(z , s(t, z )) — t. The control U constructed 
corresponding to this choice of s is the desired control. It is clearly in L 2 ([0,t];M.) 
since it is uniformly bounded. 
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We now return to the case of = where = 0. Let (X^,Y^) be the 

solution with control Ut = —2XtYt — 1 and initial condition (Xq ,Y^) = z*. Since 
this is in fact a flow, its solutions also exist backward in time. It is clear from the 
choice U that =1^0. For future reference, define z\ — {X u _ l , Y^). Hence by 
the first part of the proof, there exists a T* so that for any s > T* there exists a 
control U so that (Xg, Ytf) = z and (X^,Y S U ) = z±. Thus if we define 

- = (u t ifte[o,«] 

* \-2X t Y t -l ifte(s,s + l] ' 

By the choice of Z\, (X^ , Yf) — z* as required if we set T = s + 1. 

Lastly we observe that none of the control paths employed arc identically zero 
for all time. Hence Lemma fl 1 .21 implies that Mt, Zq (U) is non-degenerate. □ 

Combining the results from this section, we get the following Proposition which 
is the analogue of Proposition lll.il 

Proposition 11.6. If cr y > but a x = 0, then for any z* = (a;*,j/*) <E R 2 with 
< there exists a so that for any t > T* one has pt{z, z*) > for any z £ M 2 . 
Furthermore, if zq = (xo,yo) £ ^ 2 with xq < and Z\ — {x\,yi) £ R 2 with x\ > 
than one has Pt{zo, Z\) = for any t > 0. 

Proof of Proposition ] ll.lft As already outlined, the positivity claim follows from 
Theorem [4J because Proposition 111.51 guarantees the existence of a control with 
the needed properties. The fact that pt(zQ, Z\) = when x\ is strictly positive also 
follows from Theorem 21 provided we can show that there is no control which moves 
one from zq to Z\. To see this, observe that except for the fixed point at the origin, 
the vector field for any control always points toward the left half-plane along the y- 
axis. Hence it is impossible to leave the left half-plane. The fact that pt{zo, z-y) = 
when x\ = follows from the strict positivity of p t (zQ,zi) when X\ < and from 
the continuity of pt(zo, z±). □ 

11.3. Positivity of the Invariant Measure: Proof of Theorem [3l Assuming 
that a y > 0, we know that Pt(z, •) is absolutely continuous with respect to Lebesgue 
measure and has a smooth density. Hence if /j, is an invariant measure (and therefore 
we have /i = fiPt for any t > 0), we see that fi also has a smooth density, m, with 
respect to Lebesgue measure. The invariance implies that for all z £ R 2 and t > 

(50) m(z) = / p t {z',z)m{z')dz' 

Jr 2 

where pt is the density of P t . 

Let z be a point such that for some t > 0, p t (z',z) > for all z' £ M 2 . Since 
m integrates to one and is smooth, there must exist some open set A such that 
m(z') > for all z' £ A. Combining this observation with ([50]), we get 

m(z)> / pt{z', z)m(z')dz' > . 
J A 

Hence we deduce that m(z) is positive at any point z which satisfies the stated as- 
sumption. Applying this result to the information on the positivity of p t in Propo- 
sition 111.11 and Proposition lll.6[ we obtain the conclusions about the positivity of 
m(z) in Theorem 03 
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To deduce the statements that 771(2) = for z = (x,y) with x > if o~ y > and 
a x = 0, we use Proposition 1 1 1 . 6[ which states that if w £ H+ = {z = (x,y) £ M 2 : 
x > 0}, then pt((, w) — for all t > if £ £ where H- is defined to be the 
complement of H + . This implies that if z £ H + , then 

m(z) — / Pt{w , z)m{w)dw . 
Jh + 

Integrating this expression over H + and interchanging the order of integration, we 
get 

/ m(z)dz = / (f>{w)m(w)dw . 
Jh+ jh+ 

where 4>{w) = f H+ pt(w,z)dz. Since m{z),(f>{z) > for all z, this implies that for 
Lebesgue-almost-every z £ H+, either m(z) = or 4>{z) = 1. Yet from Propo- 
sition 111.61 we know that given any w £ H- , there exists a time t > so that 
Pt{z, w) > for all z £ M. 2 . Because z t-t pt(z, w) is continuous, we know (j>(w) < 1 
for every w £ H + . This implies that m(z) = for almost every z £ H + . Since 
m(z) is continuous, this forces m(z) — for all z £ H + . This completes the proof 
of Theorem [3l 



12. MlNORIZATION AND GEOMETRIC ERGODICITY 

We now establish the minorization condition we need to complete the proof of 
Theorem [2 Specifically, we seek a probability measure v and positive constants a, 
R, and T so that 

(51) inf P T (z, ■ ) > av( ■ ) 

{z£R 2 :\z\<R} 

and R > Kt where Kt is the constant from Lemma 16.11 This condition is a 
localized version of the classical Doeblin condition and central to the theory of 
Harris chains |Har561 IMT931 IHM08] . While the Lyapunov condition ensures the 
existence of an invariant measure and guarantees sufficiently rapid returns to the 
"center of phase space" to produce geometric mixing, the minorization condition 
ensures the existence of probabilistic mixing. 

To summarize our current situation, we pause to prove the following intermediate 
result. 



Theorem 5. If the minorization condition holds from (|5ip . then the Markov semi- 
group Pt generated by ([T]) satisfies the conclusions of Theorem^ 

Proof of Theorem^ By Theorem 1.3 in |HM08j . there exist constants a £ (0, 1) 
and /3 > such that p ^{^iPt , ^Pt) < o^P/si^ii ^2)- Results such as this are quite 
classical. Other proofs can be found, for example, in |MT93j . Combining this 
estimate with Proposition 16.21 immediately implies that for any n £ {0} U N 

(52) p^iP nT ,^PnT) < a^ppifuPr,^) < a n ( 1 + ^_ Kt ) Po {p, u p 2 ) . 
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To extend this estimate to an arbitrary t > 0, wc define a nonncgative integer n 
and r G (0,1) so that t = nT + t and observe that 

_ /1 + /3Kt\ 

Pf3{(llP t ,(l2Pt) = PpifJ-lPrPnT , PzPrPnT) < a ^ J p (Ml P T , PlPr) 

<ot y Jpo{Jii,fi2) < a T y ^ )po(pi,P2) 

As noted in Rcniark l6.3l for any j3' > there exists a constant C so that pp (vi, 1/2) < 
Cpp{v\, V2) for all probability measure i^. This completes the proof. □ 

12.1. Minorization when a x > 0. Now since for each t > 0, (z,z') H> p t (z,z') 
is continuous and everywhere positive, it is elementary that for any R > there 
exists a positive constant a = a(R, t) > so that inf{p t (z, z') : z, z' <E |z|, |z'| < 
R} > a. The minorization condition follows immediately from this, since for any 
measurable set A 

P t (z,A)= / p t {z,z')dz' > aLeb{An B R (0)) = aLeb(B R (0))is(A) 

J A 

where Leb is Lebesgue measure and v(A) = Leb(A n B R (0)) /Lcb(B R (0)). 

12.2. Minorization when a x = 0. We now state and prove a lemma which shows 
that the needed minorization condition follows quickly from continuity and a re- 
laxed positivity assumption. In Section fll.21 this relaxed positivity assumption was 
shown to hold by using a very explicit control theory argument coupled with some 
stochastic analysis. 

Lemma 12.1. Let P(z, dz') be a Markov transition kernel on M. d such that P(z, dz') = 
p(z, z')dz' with p: R d x R d ->• [0,oo) jointly continuous. If there exists z* G M and 
a compact set B such that for all z G B, p(z, z*) > there exist a G (0, 1) and a 
probability measure v such that 

inf P(z, - )>a . 

Proof. Since z t— > p{z,z*) is continuous, it achieves its minimum on the compact 
set B. Since p(z,z*) is strictly positive for all z G B, we know that for some 
a G (1/2,0), p(z, z*) > 2a for all z G B. By the joint continuity of p, for each z G B 
there exists a S z > be such that p(y,x) > a for all (y,x) G B$ z {z) x B^z*). 
Since {Bg z (z): z G Z?} is an open cover of the compact set B, we can extract a finite 
subcover. Let K be the collection of points z associated with this finite subcover. 
If 5 = xnva{8 z : z G K}, then S > since 8 Z > and K is finite. Now since 



BC (J B ff .(i 



we have that p(y, x) > a for all (y, x) G BxBs(z*). Then P(y, A) > aLeb(Bs^ z ){z))v{A) 1 
where v(A) — Leb(An Bg(z*)) /Leb(Bg/ z \(z)) and Leb is Lebesgue measure on R d , 
since 



P(y,A)= / p(y,x)dx> / p(y, x)dx > aLeb(A n B s (z*)) . 

JA J Ac\B s (z*) 



□ 
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13. Conclusion 



We describe a general methodology for building a Lyapunov function in a setting 
where the global stability of the systems requires flux of probability into regions 
which are clearly dissipative from the rest of phase space. We are most interested 
in problems, like the example considered here, where the noise plays an essential 
role in creating this transport is some regions. The algorithm makes use of local 
Lyapunov functions, which are constructed as solutions to Poisson equations in 
different regions, and are then patched together to form one global Lyapunov func- 
tion. We apply these techniques to one specific example in the plane to illustrate 
how the addition of noise gives rise to an invariant probability measure for a system 
whose purely deterministic dynamics exhibit instability. Furthermore, our resulting 
"super" Lyapunov function enables us to extract a stronger convergence than what 
is usually proved in the Harris-chain setting — indeed, an exponential convergence 
independent of initial condition — to this equilibrium measure. En route to proving 
this convergence, we employ explicit control-theoretic constructions and we rely on 
the tools of Malliavin calculus. It is our hope that the simple and specific appli- 
cations of control theory and Malliavin calculus in our model problem will be of 
independent interest. 

Of course, further work remains to be done: the application of these methods 
to other examples, for instance, and the development of more general theorems 
about noise-induced stabilization. In particular, in the planar system we consider, 
patching the local Lyapunov functions turns out to be one of the most delicate and 
important parts of the proof. Therefore, it would be especially interesting to find 
more general approaches to the problem of patching Lyapunov functions, and more 
general conditions under which it can be done successfully. 

When our construction works, it is likely to produce a Lyapunov function which 
provides strong control over the excursions towards infinity and a nearly sharp rate 
for the convergence to equilibrium. However, the construction of such a function 
is laborious. It would be interesting to obtain a simpler "partial fluid limit" which 
captures only the minimal stochasticity at infinity needed to stabilize the system. 
This may arise as an extension to our work in the direction of |DW94[ IHKM02 , 
FMMP08, Mey08 . Such an approach might allow simpler proofs of stabilization 
without necessarily proving the existence of strong Lyapunov function. 



Proposition A.l. Suppose f £ C(K) is a non-increasing function and that <p(t) 
and ip(t) are C 1 functions on R satisfying (f)(0) = ip(0) and 4>'(t) < f(<p(t)), ip'(t) — 
f(ip(t)) for allt>0 then <f>(t) < if)(t) for all t > 0. 

Proof of Proposition \A. l[ For all < r < t, we have that 



Let T\ = inf{i > : ijj(t) — <p(t) < 0}. Suppose for contradiction that T\ < oo. 
Then by continuity, ip(Ti) — 4>{T\) = and there exists T 2 £ (Ti,oo) such that for 
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and 




which implies 
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all t £ (Ti,T 2 ), ip(t) - 4>(t) < 0. Then for all t £ (Ti,T 2 ), 

(53) - 0(*) > (V-(Ti) - HTi)) + f (f(m) ~ f(<f>(s)))ds . 

Now since / is non-increasing, ip(t) < <j>(i) implies that f(ip(t)) > f(<p(t). Hence 
this combined with ((53|) implies that for all t £ (Ti,T 2 ), tp(t) - <p(t) > 0. This is a 
contradiction. Hence T\ must be infinite and 4>(t) < ip(t) for all t > 0. □ 
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